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This application is based upon an earlier filed U.K Patent Application 
5 No. 9301268.0, filed January 22, 1993, which, in turn, is a continuation-in-part of 
U.K. Patent Application No. 9216383.1, filed July 31, 1992, which, in turn, is a 
continuation-in-part of U.K. Patent Application No. 9210810.9, filed May 21, 1992, 
which, in turn, is a continuation-in-part of U.K. Patent Application No. 9209648.6, 
filed May 5, 1992, which, in turn, is a continuation-in-part of U.K. Patent Application 
10 No. 9207013.5, filed March 30, 1992, which, in turn, is a continuation-in-part of U.K. 
Patent Application No. 9205541.7, filed March 13, 1992, which, in turn, is a 
continuation-in-part of parent U.K. Patent Application No. 9205058.2, filed March 9, 
1992, the benefit of the filing dates of which is hereby claimed pursuant to 

35U.S.C. § 119. 
15 Field of the Invention 

The present invention relates generally to the field of imaging and, more 
particularly, to the imaging of nerve tissue and other diffusionally anisotropic 
structures. 

Background of the Invention 
20 Although many techniques have been developed for locating and viewing the 

brain, spinal cord, and spinal roots within the spinal canal, hitherto there has not been 
a successfiil method for viewing the peripheral, autonomic, and cranial nerves. These 
nerves, collectively referred to herein as peripheral nerves, commonly travel through 
and along bone, muscle, lymphatics, tendons, ligaments, intermuscular septa. 
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collections of fatty tissues, air and fluid spaces, veins, arteries, joints, skin, mucous 
membranes and other tissues. The relatively small size of peripheral nerves, as well as 
their close proximity to other tissue of comparable size and shape, makes them 

difficult to locate and identify. 
5 The examination of peripheral nerves is further complicated by the complexity 

of many such neural structures, such as the brachial plexus, lumbar plexus, and sacral 
plexus. These structures include bundles of nerves that may join together, separate, 
rejoin, intermix, and resegregate, forming intricate three dimensional patterns. A 
compression or irritation of a small area of nerve within such a plexus (e.g. in the 

10 shoulder) can cause pain, numbness, weakness or paralysis at some distant site (e.g. in 
one finger). Even when a surgeon attempts to expose the plexus for direct inspection, 
the anatomic complexity can prove overwhelming, rendering diagnosis inconclusive 
and surgery difficult and dangerous. 

Radiologic methods employing, for example, X-rays, have been developed to 

15 generate tissue specific images of various physiological structures including bone, 
blood vessels, lymphatics, the gastrointestinal tract, and the tissues of the central 
nervous system; Due in part to the neural characteristics noted above, however, these 
techniques have not been successfully used to generate suitable clinical images of 
peripheral nerves. 

20 Typically, the position of peripheral nerves in radiologic images has been 

inferred by reference to more conspicuous, non-neural structures such as tendons, 
vessels, or bone. For example, by producing an X-ray image of a region of the body 
through which a nerve of interest passes, non-neural structures can often be readily 
identified. Then, the locations of peripheral nerves in the region can be inferred from 

25 standard reference information about human anatomy. Due to the variability of nerve 
position from one individual to another, however, this technique is of limited value. 

One approach of particular interest tfiat has been used to image physiological 
structures is magnetic resonance imaging (MRI). By way of introduction, MRI 
involves the exposure of tissue to a variety of different magnetic and radio-frequency 

30 (rf) electromagnetic fields. The response of the specimen's atomic nuclei to the fields 
is then processed to produce an image of the specimen. 

More particularly, the specimen is initially exposed to a polarizing magnetic 
field. In the presence of this field, nuclei exhibiting magnetic moments (hereinafter 
referred to as spins) will seek to align themselves with the field. The nuclei precess 

35 about the polarizing field at an angular frequency (hereinafter referred to as the 
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Larmor frequency) whose magnitude depends upon both the field's strength and the 
magnetogyric constant of the specific nuclear species involved. 

Although the magnetic components of the spins cancel each other in a plane 
perpendicular to the polarizing field, the spins exhibit a net magnetic moment in the 
5 direction of the polarizing field. By applying an excitation field perpendicular to the 
polarizing field and at a frequency near the Larmor frequency, the net magnetic 
moment can be tilted. The tilted magnetic moment includes a transverse component, 
in the plane perpendicular to the polarizing field, rotating at the Larmor frequency. 
The extent to which the magnetic moment is tilted and, hence, the magnitude of the 
10 net transverse magnetic moment, depends upon the magnitude and duration of the 
excitation field. 

An external return coil is used to sense the field associated with the transverse 
magnetic moment, once the excitation field is removed. The return coil, thus, 
produces a sinusoidal output, whose frequency is the Larmor frequency and whose 
15 amplitude is proportional to that of the transverse magnetic moment. With the 
excitation field removed, the net magnetic moment gradually reorients itself with the 
polarizing field. As a result, the amplitude of the return coil output decays 

exponentially with time. 

Two factors influencing the rate of decay are known as the spin-lattice 

20 relaxation coefficient T^ and the spin-spin relaxation coefficient T2. The spin-spin 
relaxation coefficient T2 represents the influence that interactions between spins have 
on decay, while the spin-lattice relaxation coefficient T^ represents the influence that 
interactions between spins and fixed components have on decay. Thus, the rate at 
which the return coil output decays is dependent upon, and indicative of, the 

25 composition of the specimen. 

By employing an excitation field that has a narrow frequency band, only a 
relatively narrow band within a nuclear spicies will be excited. As a result, the 
transverse magnetic component and, hence, return coil output, will exhibit a relatively 
narrow frequency band indicative of that band of the nuclear species. On the other 

30 hand, if the excitation field has a broad frequency band, the return coil output may 
include components associated with the transverse magnetic components of a greater 
variety of frequencies. A Fourier analysis of the output allows the different 
fi-equencies, which can be indicative of different chemical or biological environments, 
to be distinguished. 

35 In the arrangement described above, the contribution of particular spins to the 

return coil output is not dependent upon their location within the specimen. As a 
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result, while the frequency and decay of the output can be used to identify 
components of the specimen, the output does not indicate the location of components 
in the specimen. 

To produce such a spatial image of the specimen, gradients are established in 
5 the polarizing field. The direction of the polarizing field remains the same, but its 
strength varies along the x, y, and z axes oriented with respect to the specimen. By 
varying the strength of the polarizing field linearly along the x-axis, the Larmor 
frequency of a particular nuclear species will also vary linearly as a Sanction of its 
position along the x-axis. Similarly, with magnetic field gradients established along 

10 the y-axis and z-axis, the Larmor frequency of a particular species will vary linearly as 
a fiinction of its position along these axes. 

As noted above, by performing a Fourier analysis of the return coil's output, 
the frequency components of the output can be separated. With a narrow band 
excitation field applied to excite a select nuclear species, the position of a spin relative 

15 to the xyz coordinate system can then be determined by assessing the difference 
between the coil output frequency and the Larmor frequency for that species. Thus, 
the MRI system can be constructed to analyze frequency at a given point in time to 
determine the location of spins relative to the magnetic field gradients and to analyze 
the decay in frequency to determine the composition of the specimen at a particular 

20 point. 

The generation and sensing of the fields required for proper operation of an 
MRI system is achieved in response to the sequential operation of, for example, one 
or more main polarizing field coils, polarizing gradient field coils, rf excitation field 
coils, and return field coils. Commonly, the same coil arrangement is used to generate 

25 the excitation field and sense the return field. A variety of different sequences have 
been developed to tailor specific aspects of MRI system operation, as described, for 
example, in U.S. Patent No. 4,843,322 H:Glover); U.S. Patent No. 4,868,501 
(ConoUy); and U.S. Patent No. 4,901,020 (Ladebeck et al.). 

One application of conventional MRI systems is in the production of 

30 angiograms, or blood vessel images. Various different pulse sequences and 
processing techniques have been developed for use in MRI angiography, as described 
in, for example, U.S. Patent No. 4,516,582 (Redington); U.S. Patent No. 4,528,985 
(Macovski); U.S. Patent No. 4,647,857 (Taber); U.S. Patent No. 4,714,081 
(Dumoulin et al.); U.S. Patent No. 4,777,957 (Wehrli et al.); and U.S. Patent No. 

35 4,836,209 (Nishimura). 
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As will be appreciated, blood vessels are readily differentiated from 
surrounding tissue by the pulsatile flow of blood therethrough. MRI angiography 
exploits this distinguishing characteristic to generate images of the blood vessels in 
various ways. For example, if the excitation field is pulsed at systole and diastole, the 
5 contribution of blood flow to the return field will differ, while the contribution of 
static tissue and bone to the return field will be the same. By subtracting one return 
fi-om the other, the static component cancels, leaving only the contribution fi-om the 
blood vessel. 

Unfortunately, because peripheral nerve does not exhibit the flow- 

10 distinctiveness of blood vessels, MRI angiography systems and pulse sequences can 
not be used to generate suitable images of peripheral nerve. Further, conventional 
MRI systems and sequences used for general imaging of tissue and bone do not 
provide acceptable results. Given the poor signal-to-noise (S/N) ratio of the return 
signals (e.g., on the order of Ix to 1.5x) and the small size of the nerve, the 

15 conspicuity of imaged nerves relative to other tissue is collectively rendered so poor 
as to be diagnostically useless. 

One technique proposed for use in enhancing the utility of MRI systems in 
imaging neural tissue involves the use of pharmaceutical agents to enhance the 
contrast of neural tissue relative to surrounding tissue in the images produced. As 

20 described in PCX Patent^pplic^ion No. PCX EP 91/01780 (Filler et al., WO 
92/04916), published on^; ^^cl i 1 5, 1992, a two-part contrast agent, such as wheat 
germ agglutinin or dextrin-magnetite, is injected so that it is subsequently taken up, 
and transported, by the nerve of interest. Xhe first part of the agent promotes neural 
uptake, while the second part of the agent has the desired "imageable" property. 

25 Xhe agent is injected into muscle and undergoes axoplasmic flow in the nerve 

supplying that muscle, tagging the nerve in subsequently generated images of the 
specimen If MRI is used, the second part of the agent is selected to have a 
magnetically active (e.g., ferrite) component. An agent having a high nuclear density 
can, however, be used to increase the contrast of the nerve upon X-ray or computed 

30 tomography (CX) examination, while a radioactive (e.g. positron emitting) element 
can be used to enhance visibility during positron emission tomography (PEX) 
scanning. 

Xo illustrate the effectiveness of contrast agents in imaging nerve, reference is 
had to FIGURES 1-5. In that regard, FIGURE 1 is a diagram of a transverse section 
35 of the upper forearm FA of a rabbit. Xhe forearm includes the triceps muscle XM, 
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ulnar nerve UN, brachial veins BV, median nerve MN, radial nerve RN, humerus H, 
cephalic vein CV, and biceps muscle BM. 

FIGURES 2A and 2B illustrate spin-echo MR images of such a section, using 
a ferrite contrast agent, produced by a conventional MRI system at six-hour intervals. 
5 Although some of the larger structural elements are readily identified, the location of 
some objects appears skewed. More particularly, the humerus marrow appears 
shifted relative to the humerus H, as do ligaments L, and fat F between the biceps or 
triceps. In addition, smaller neural structures are difficult to distingu^. 
^ Several approaches are available, however, to attempt te^t d e ntity -rter^es in the 

10 images generated. For example, as shown in FIGURE 3, if a short tau inversion 
recovery (STIR) sequence of the type described in Atlas et al., STIR MR Imaging of 
the Orbit, 151 AM. J. ROENTGEN. 1025-1030 (1988) is used, the humerus marrow 
disappears from the image as does, more importantly, certain ambiguous, apparently 
non-neural structures adjacent the median nerve MN. Thus, as shown in the enlarged 

1 5 image of the region including the median nerve MN and ulnar nerve UN, provided in 
FIGURE 4, the median nerve MN is visible. 

Similarly, even when the contrast agent images of FIGURES 2A and 2B are 
enlarged to better illustrate the region including the median nerve MN, as shown in 
FIGURES 5A and 5B, respectively, the nerves are distinguishable to a highly skilled 

20 observer. More particularly, transport of the ferrite contrast agent during the six-hour 
interval between the generation of images 4A and 4B results in a loss of intensity in 
the MN relative to the non-neural structure adjacent median nerve MN, Given this 
observation and the STIR-based assessment, the median nerve MN can, thus, be 
identified. 

25 The use of contrast agents, while promising, does have certain limitations. For 

example, there is an increasing preference to avoid the use of invasive technologies in 
medicine whenever possible. Further, contrast agents generally can be used to image 
only a single nerve or nerve group. Of perhaps greatest importance, the contrast 
agents employed typically reduce the intensity of the imaged nerve. Since nerves are 

30 already difficult to see in current MRI images, the impact of the contrast agent upon 
the image can be diflFicult to interpret, as illustrated by the discussion of FIGURES 2- 
5 above. 

In another application, MRI has been used, without contrast agents, to map 
non-peripheral, white matter nerve tracts in the brain. The white matter tracts extend 
35 through gray matter tissue in the brain and exhibit relatively high anisotropic diffusion. 
More particularly, given their physical structure (i.e., axonal pathways surrounded by 
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myelin sheaths), water mobility along the white matter tracts is relatively high, while 
water mobility perpendicular to the tracts is low. The surrounding gray matter does 
not, however, exhibit this same anisotropy. 

A technique for MRI-based mapping of white matter nerve tracts that exploits 
5 this characteristic of neural tissue is described in Douek et al.. Myelin Fiber 
Orientation Color Mapping, BOOK OF ABSTRACTS, SOCIETY OF MAGNETIC 
RESONANCE IN MEDICINE, p. 919 (1991). Basically, in addition to the fields and 
gradients described above, this process involves the use of a pair of field gradient 
pulses (hereinafter referred to as diffusion gradients), oriented perpendicular and 

10 parallel to the white matter tracts to be imaged. The effect of a pulsed gradient is to 
change the phase of the received signal from all of the spins. For stationary spins the 
effect of the two diffusion gradients cancels out. In contrast, spins moving from one 
spatial position to another in the time between the two diffusion gradients experience 
changes in the frequency and phase of the spin magnetization with the net effect being 

15 a reduction in the received signal. The signal reduction is greatest for those spins that 
diffuse the greatest distance between the two pulsed gradients. 

As noted above, given the anisotropic nature of the tracts, water will diffuse 
freely along a tract, but is restricted in it motion perpendicular to the tract. When the 
diffusion gradient is aligned with the tract there is thus a greater reduction in signal 

20 than when the diffusion gradient is aligned perpendicular to the tract. Because this 
phenomenon is not exhibited by the surrounding gray matter tissue, the white matter 
tracts can be identified. 

Anisotropic diffusion is also a recognized characteristic of peripheral nerve, as 
indicated in Moseley et al., Anisotropy in Diffusion-Weighted MRI, 19 MAGNETIC 

25 RESONANCE ON MEDICINE 321 (1991). The Douek et al. technology, however, 
does not distinguish peripheral nerve from muscle and other tissue for a number of 
previously unrecognized reasons. First, wfiile the size and structure of the white 
matter tracts ensure that the resultant signals will be sufficiently strong for imaging, 
peripheral nerve is considerably smaller and more difficult to distinguish. Second, 

30 unlike the white matter tracts, peripheral nerve is commonly surrounded by muscle 
and fat, both of which impair the ability of the Douek et al. system to image nerve. 

By way of elaboration, given its fibrous structure, muscle also exhibits 
diffiisional anisotropy, as recognized in Moseley et aL, Acute Effects of Exercise on 
Echo-Planar T2 and Diffusion-Weighted MRI of Skeletal Muscle in Volunteers, 

35 BOOK OF ABSTRACTS, SOCIETY OF MAGNETIC RESONANCE IN 
MEDICINE 108 (1991). As a result, the simple anisotropic analysis of Douek et al. is 
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unable to distinguish peripheral nerve and muscle. While fat is isotropic and, 
therefore, distinguishable from nerve, it also impairs the imaging of peripheral nerves. 
Specifically, the relative signal strength of fat returns to neural returns is so high as to 
render peripheral nerves unidentifiable in images produced. 
5 As will be appreciated from the preceding remarks, it would be desirable to 

develop a method for rapidly and non-invasively imaging a single peripheral nerve, or 
an entire neural network, without resort to contrast agents. The images generated 
should be sufficiently detailed and accurate to allow the location and condition of 
individual peripheral nerves to be assessed. It would fiarther be desirable to provide a 
10 system that processes neural images to enhance the information content of the images, 
diagnose neural trauma and disorders, and inform and control the administration of 
treatments and therapy. 

Summary of the Invention 
The present disclosure relates to a new method, which quite remarkably, is 

15 capable of generating a three dimensional image of an individual patient's nerves and 
nerve plexuses. The image can be acquired non-invasively and rapidly by a magnetic 
resonance scanner. These images are acquired in such a way that some embodiments 
of the invention are able to make all other structures in the body including bone, fat, 
skin, muscle, blood, and connective tissues tend to disappear so that only the nerve 

20 tree remains to be seen. A plurality of the nerves passing through a given imaged 
region may be observed simultaneously, thus alleviating any ambiguity of nerve 
identification which might arise were only a single nerve imaged as with some contrast 
agent techniques. 

The invention is based on the discovery of a method of collecting a data set of 
25 signal intensities with spatial coordinates which describes the positions of the nerves 
within any two dimensional cross section of a living mammal or within any three 
dimensional data acquisition space. There exist a large number of pulse sequences 
capable of controlling or operating a magnetic resonance imaging apparatus and each 
of which accomplishes some preferred image optimization. Previously, however, no 
30 simple (single) or complex (double or multiple) pulse sequence has been able to 
increase the relative signal intensity of nerve so that it is brighter than all other tissues 
in the body or limb cross section. Surprisingly, the inventors have discovered that 
there are certain novel ways of assembling complex pulse sequences, wherein even 
though the simple components of the sequence decrease the signal-to-noise ratio of 
35 nerve or decrease the signal strength of nerve relative to other tissues, the fiilly 
assembled complex sequence actually results in the nerve signal being more intense 
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than any other tissue. In this fashion, the image conspicuity of nerve is greatly 
increased. 

Thus, a first aspect of the present invention provides a method of selectively 
imaging neural tissue of a subject without requiring use of intraneural contrast agents, 
5 the method comprising subjecting part of the subject anatomy to magnetic resonance 
imaging fields, detecting magnetic resonance and producing an image of neural tissue 
fi-om said detected resonance so that a nerve, root, or neural tract of interest in said 
image can be visually differentiated from surrounding structures. 

A second aspect of the present invention provides a method of selectively 

10 imaging neural tissue of a subject, the method comprising subjecting part of the 
subject anatomy to magnetic resonance imaging fields adapted to discriminate 
anisotropy of water diffusion or other special characteristic of neural tissue, detecting 
magnetic resonance to produce an electronic signal in accordance with said resonance 
and producing an image of neural tissue from said electronic signal. 

15 The invention also provides an apparatus for selectively imaging neural tissue 

of a subject without requiring the use of neural contrast agents, the apparatus 
comprising means for subjecting part of the subject anatomy to magnetic resonance 
fields, means for detecting magnetic resonance to produce an electronic signal in 
accordance with said resonance, and means for producing an image of neural tissue 

20 from said electronic signal so that a nerve, root, or neural tract of interest in said 
image can be visually differentiated from surrounding structures. 

The invention also finds expression as an apparatus for imaging neural tissue 
of a subject, the apparatus comprising means for subjecting part of the subject 
anatomy to magnetic resonance fields adapted to discriminate anisotropy of water 

25 diffusion, means for detecting magnetic resonance to produce an electronic signal in 
accordance with said resonance and means for producing a selective image of neural 
tissue of interest from said electronic signal. ^ 

Brief Description of the Drawings 
30 The foregoing aspects and many of the attendant advantages of this invention 

will become more readily appreciated as the same becomes better understood by 
reference to the following detailed description, when taken in conjunction with the 
accompanying drawings, wherein: 

FIGURE 1 is a diagram of a transverse section of the upper forearm of a 
35 rabbit illustrating various neural and non-neural structures; 
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FIGURES 2A and 2B are images of the upper forearm of a rabbit, of the type 
depicted in FIGURE 1, produced using an MRI system at two spaced-apart times 
after the forearm was injected with a ferrite contrast agent; 

FIGURE 3 is another image of the upper forearm of a rabbit produced using 
5 an MRI system employing a short tau inversion recovery (STIR) spin-echo sequence; 

FIGURE 4 is an enlargement of a portion of the image of FIGURE 3 
associated with a peripheral nerve of interest; 

FIGURES 5A and 5B are enlargements of a portion of the images of 
FIGURES 2 A and 2B, respectively, associated with a peripheral nerve of interest; 
10 FIGURE 6 is a block diagram of a neurography system, constructed in 

accordance with this invention, coupled to a plurality of other systems designed to 
provide information to the neurography system and to implement, for example, neural 
diagnoses, therapy, surgery, and training; 

FIGURE 7 is a functional chart of the operation of the neurography system of 
15 FIGURE 6; 

FIGURE 8 is an illustration of the various components included in the 
neurography system of FIGURE 6; 

FIGURES 9 and 10 are flow charts depicting one way in which the 
neurography system of FIGURE 8 may be used to generate neurograms; 
20 FIGURES llA through IIF illustrate one sequence of pulses suitable for use 

in producing diagnostically suitable images from the neurography system of 
FIGURE 6; 

FIGURE 12 is another image of the upper forearm of a rabbit produced by an 
embodiment of the neurography system employing fat suppression; 
25 FIGURES 13 A and 13 B are additional images of the upper forearm of a 

rabbit produced by an embodiment of the neurography system employing gradients 
perpendicular and parallel, respectively, to the anisotropic axis of nerve being imaged; 

FIGURES 14A through 14D are images of the upper forearm of a rabbit 
produced employing gradients of 0, 3, 5, and 7 Gauss/centimeter, respectively; 
30 FIGURES ISA through 15C are images produceable by the neurography 

system with zero, perpendicular, and parallel gradients, while FIGURE 15D is an 
image based upon the images of FIGURES 15B and 15C, referred to herein as a 
subtraction neurogram 

FIGURE 16 vector length image of the brain produced using the neurography 
35 system of FIGURE 8; 
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FIGURE 17 is an arctan image of the brain produced using the neurography 

system of FIGURE 8; 

FIGURES 18A through 18D are images of a rabbit forearm produced using 
the neurography system of FIGURE 8, and illustrating the influence of the TE 
5 sequence upon the images produced; 

FIGURE 19 illustrates a splint employed in the neurography and medical 
systems of the present invention; 

FIGURES 20, 21, and 22 are illustrations of images of a human sciatic nerve 
produced using the neurography system of FIGURE 8, with FIGURES 20 and 21 
10 illustrating the ability of the system to image nerve fascicles (in two cross-sectional 
scales) and FIGURE 22 illustrating an axial projection of the nerve; 

FIGURE 23 is a diagram of a cross-section of a vertebra, illustrating the types 
of structure present in one neurography application; and 

FIGURE 24 is a schematic illustration of a surgical system constructed in 
15 accordance with this invention for use with the neurography system of FIGURE 8. 

Detailed Description of the Preferred Embodiment 
Referring now to FIGURE 6, a neurography system 10 is shown as one 
component of a broader medical system 12. Unlike prior art arrangements, system 10 
quickly and non-invasively generates accurate images showing the pattern of 
20 individual peripheral nerves, or entire nerve trees, without the use of contrast agents. 
The system is designed to allow such images, hereinafter referred to as neurograms, 
to be displayed in two-dimensions, illustrating neural cross sections in the specimen 
under examination, or in three-dimensions. The images may selectively exclude all 
other structures within the specimen, or may illustrate the physical relationship of 
25 other structures relative to the nerves for reference. 

1 . Medical Svstem Overview 
As shown in FIGURE 6, the neurography system 10 included in medical 
system 12 includes four basic components: MRI system 14, processing system 16, 
input system 18, and output/display system 20. In the preferred arrangement, the 
30 MRI system 14 is a conventional MRI system modified for use in collecting image 
data of a patient P under examination. The processing system 16 responds to 
operator inputs applied via input system 18 to control MRI system 14 and process its 
output to display the resultant neurograms at system 20. As will be described in 
greater detail below, system 16 employs a variety of different imaging protocols, 
35 alone or in combination, to ensure that the images produced are of a quality 
heretofore unachieved. 
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The medical system 12 includes a number of components that supplement the 
imaging information produced by system 10 and/or use that information for a variety 
of purposes. For example, an auxiliary data collection system 22 may be included to 
collect image information about non-neural structures, such as blood vessels and 
5 bone, in the imaged region of patient P. This information can then be used to 
suppress and/or enhance the appearance of those structures in the neurograms 
produced by system 10. 

A diagnostic system 24, included in system 12, may be used to analyze the 
images produced by system 10. Given the high resolution, detail, and accuracy of 

10 neurograms produced by system 10, system 24 can be programmed to analyze neural 
pathway information to detect discontinuities associated with, for example, neural 
compressions, injuries, and tumors. System 24 provides outputs indicative of the 
location of discontinuities and may, by consultation with a database of image 
information associated with clinically assessed abnormalities, provide an indication of 

15 the nature and magnitude of an imaged discontinuity. These outputs can be used for 
diagnosis, or applied as feedback to system 10 to refine a region of interest (ROI) 
under examination in patient P. 

Medical system 12 may also include a therapeutic system 26 and surgical 
system 28. Systems 26 and 28 employ information about the patient's neural structure 

20 from system 10 to assist in the proper administration of a desired therapeutic or 
surgical operation. For example, the information may be used to guide a robotic 
stylus to a damaged neural site for treatment or to allow an operation on non-neural 
structure to be performed without damage to the patient's peripheral nerves. The 
systems 26 and 28 may operate independent of physician control or may simply 

25 provide the physician with real-time feedback concerning the relationship between an 
operation being performed and the patient's neural structures. 

A training and development system 36 is included in the medical system 12 for 
a variety of different purposes. For example, the training system 30 may be used to 
demonstrate the anatomy of various neural structures, along with their positional 

30 relationship to non-neural patient structures. This information has great educational 
value given the extremely limited ability of prior art techniques, including direct 
examination, to provide detailed anatomical information. Training system 30 may also 
be designed to analyze the effectiveness of neurography system 10 and provide 
feedback used to control the pulse sequences and other operational parameters of 

35 system 10. 


UOFW\8060AP.DOC 



-13- 


0^ 


As one final component, medical system 12 may include a host processing 
system 32 in addition to, or in place of, separate processing systems in the other 
components of system 12. Although not separately shown in FIGURE 6, system 32 
includes a central processing unit (CPU) coupled to the remainder of system 12 by 
5 input/output circuits. Memory is provided to store software instructions, used to 
control the operation of the CPU and, hence, the various components of system 12, 
and to store image and other data collected by system 12. The use of a separate host 
processing system 32 is particularly desirable where various components of system 12 
are to be operated in interactive fashion pursuant to a single set of software 
10 instructions, 

2. The Neurography System 
Turning now to a more detailed discussion of neurography system 10, by way 
of introduction, some of the more important operational features of system 10 are 
loosely depicted in the chart of FIGURE 7. As will be described in greater detail 

15 below, system 10 may be constructed to employ one or more^these features to 
enhance the imaging ability of conventional MRI sufficiently to provide diagnostically 
and therapeutically useftil information. 

As shown, the operation of system 10 can be broken down into the broad 
steps of data collection 34, image processing and analysis 36, image display 38, and 

20 control 40. The data collection process 34 involves, for example, spin-echo 
imaging 42, which may be supplemented by one or more of the following imaging 
protocols: fat suppression 44, diffusion weighting 46, and "long T2" processing 48, 
and other protocols including magnetization transfer. Each of these protocols has 
been found to enhance the quality of images of peripheral nerve sufficiently to provide 

25 heretofore unavailable MRI neurograms. 

The data collected by process 34 is subjected to image processing and 
analysis 36, involving two-dimensional and -three-dimensional image generation 50. 
Image generation 50 may be further enhanced by miscellaneous suppression 
features 52, responsible for reducing the influence of, for example, blood vessels and 

30 patient motion, on the images produced. An image subtraction feature 54 may also be 
employed to remove all non-neural components from the images. 
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a. Neurography System Construction 
Having briefly summarized the operational aspects of neurography system 10, 
its construction and operation will now be considered in greater detail. In one 
embodiment, MRI system 14 includes an imager I of the type sold by GE Medical 
5 Systems, under the trademark SIGNA (software release 5.2). 

In that regard, as shown in FIGURE 8, the region R of the patient to be 
imaged is placed within the bore B of the MRI system imager I. As will be described 
in greater detail below, the position of region R relative to the imager may be 
stabilized by a splint 58. Splint 58 limits motion artifact, provides fiducial markers in 
10 a secondary frame of reference, and reduces the system's susceptibility to boundary 
effects that otherwise might degrade fat suppression near the boundary between skin 
and air. , > / 

MRI system 14 includes polarizing field coils 60^ responsible for exposing 
region R to the desired polarizing field. The polarizing field has a strength of, for 
15 example, 1.5 Tesla and is oriented along a z-axis. 

A tuned rf excitation coil 62 is also positioned within bore B over the region R 
under investigation. Coil 62 is provided with a pulsed rf input, in a manner described 
below, to generate the field responsible for excitation of nuclear spins in region R. 
Coil 62 is also responsible for detecting the rf return, or echo, fields generated by the 
20 spins, although separate transmit and receive coils may alternatively be used. 

The excitation coil 62 may be, for example, a solenoid or surface coil, 
configured and dimensioned to fit closely over the region R to be imaged (e.g., the 
patient's arm, leg, shoulder, chest, pelvis, head, neck or back). In a preferred 
arrangement, however, a phased array coil system is employed to increase the signal- 
25 to-noise ratio of the returns, thereby providing an improvement in the spatial 
resolution of system 14 and allowing information to be retrieved from signals that 
would otherwise have been too weak to fcwm useful images. For example, where 
peripheral nerve having a thickness on the order of 1 -2mm is to be sharply resolved, 
each array includes, for example, 4-6 individual coils, arranged in transverse and 
30 longitudinal pairs or linear paired arrays. 

Three pairs of gradient coils 64 and 66 are also positioned within the bore B 
of the imager. These coils superimpose a locational gradient of roughly one Gauss 
per centimeter upon the polarizing field over the sample region R along each of the x, 
y, and z-axes. For the sake of simplicity, however, only the z-gradient coils 64 and 66 
35 are shown in FIGURE 8. 
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In the preferred arrangement, the same coil pairs 64 and 66 are used to 
produce diffusional gradients along the desired axes, as well as the requisite locational 
gradients. Alternatively, one or more separate diffusional gradient coil pairs 68 and 
70 may be provided within the imager bore B. If the separate coil pair 68 and 70 is 
5 mounted on a movable track, substantially any desired diffusional gradient orientation 
can be achieved. The diffusional gradient is relatively strong compared to the 
locational gradients, e.g, ranging up to 10 Gauss/ centimeter or higher. 

A computer 72 and front-end circuit 74 form the processing systgm 16, input 
system 18, and output/display system 20 of neurography system shown in 
0^ 10 FIGURE 6. Computer 72 and circuit 74 cooperatively control and synchronize the 

operation of MRI system 14, as well as process and display the acquired data. 

The computer 72 is, for example, an IBM-compatible personal computer 
including a 486 processor, VGA monitor, and keyboard. An interface bus 76, 
included in circuit 74, couples computer 72 to the other components of circuit 74. 

15 A gradient pulse generator 78 included in circuit 74 produces generally 

rectangular output pulses used to establish the desired gradients in the polarizing field. 
The output of generator 78 is applied to x-, y, and z-axis gradient field amplifiers 80, 
although only the z-axis amplifier 80 is shown in FIGURE 8. As will be appreciated, 
if separate coils 68 and 70 are employed to establish the diffusional gradients, the 

20 output of generator 78 must be applied to those coils via separate amplifiers 82. 
^ Circuit 74 also includes an rf pulse generator 84, which produces rf signal^ 

pulses used in the establishment of the excitation field. In the preferred arrangement, 
the pulse generator produces an rf output suitable for use in proton MRI, although 
frequencies specific to other MRI susceptible nuclei, such as, ^ ^fluorine, carbon, 

25 ^iphosphorus, deuterium, or 23sodium, may be used. The output of generator 84 is 
amplified by a high-power rf amplifier 86 before being selectively applied to the 
excitation coil 62 by a duplexer 88. The duplexer 88 is also controlled to selectively 
steer the low level MR returns received by the excitation coil 62 to a preamplifier 90. 

A mixer 92 transforms the high frequency output of preamplifier 90 to a low 

30 frequency signal by mixing the amplified MR returns with signals from a digitally 
controlled rf oscillator 94, which also provides inputs to generator 84. The analog 
output of mixer 92 is input to a low pass filter 96 before finally being converted to a 
digital form by an analog-to-digital converter 98. The computer 72 processes the 
resultant digital inputs, which represent the response of the spins to the applied fields, 

35 to generate the desired neurograms. 
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h Neurography Sy ^tftm Operation 

Having reviewed the basic construction of the neurography system 10, its 
operation to generate the desired two- or three-dimensional neurograms will now be 
considered. To that end, FIGURES 9 and 10 depict the general sequence of steps 
performed by system 10 in the production of neurograms. These neurograms exhibit 

5 a high nerve conspicuity, which for the purpose of the ensuing discussion will be 
understood to refer to the contrast (in, for example, intensity or color) between the 
nerve and the image background. The methods described below may be used to 
produce neurographic images of substantially any region of the body, including the 
brain, for example, central nervous system (CNS) neurograms. 

10 As indicated at block 100, the operation of the system is first initialized to 

establish certain parameters of the system's operation. In that regard, the operator 
may input desired parameters via computer 72 in response to queries generated at 
start up. Because most aspects of the system's operation are controlled by software 
resident in the memory of computer 72, default initialization parameters may also be 

15 accessed. 

Although the particular parameters to be initialized may vary at the user's 
discretion, examples include the type of images to be generated (i.e., two-dimensional 
cross sections or three-dimensional projections), field of view (FOV), thickness of 
each slice imaged, pulse repetition rate (TR), number of phase encoding steps, the 
20 existence of a known axis of diffiisional anisotropy, and the strengths and orientations 
of the diffiisional gradients to be used. By way of example, the operator may select a 
two-dimensional image, a FOV of four cm by four cm, a TR of 1.5 seconds, and 256 
phase encoding steps. A discussion of anisotropic axis identification is provided 
below. 

25 Once initialization has been completed, a series of steps, corresponding to the 

data collection process 34 discussed in connexion with FIGURE 7, are performed. 
This process generally involves the control of pulse sequences used in connection with 
front end circuit 74. As will be described in greater detail below, different sets of 
pulse sequences and combinations of pulse sequences have been devised to 

30 unambiguously distinguish small peripheral nerves from neighboring structures of 
similar shape and location, including the combination of certain existing sequences 
into new groupings for use in new situations and the design of new sequences that 
incorporate optimized features for the purpose of neurographic imaging. For 
illustrative purposes, a graphic illustration of one example of a suitable pulse sequence 

35 is provided in FIGURES 1 1 A through 1 1 F. 
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i Fat Suppression 

As indicated in block 102 of FIGURE 9, a first, optional, step performed in 
the image generation process is fat suppression. Although fat represents a known 
source of interference in MRI images of bone and tissue, it was not previously 
recognized as an impediment to effective neural imaging due to the broader 
5 perception that neural MR signals were inadequate for imaging regardless of 
background composition. The value of fat suppression was discovered during the 
development of the present invention by the fortuitous use of a main field magnet 
designed for spectroscopy as part of an imaging system. 

In that regard, in MR spectroscopy, a relatively strong magnetic field is 
10 employed to increase the separation in fi-equency between signals arising fi-om 
different chemical species of the same nucleus, thereby allowing these components to 
be more easily distinguished. MRI also uses a frequency distribution (created by 
applying a field gradient) over a sample to locate spins and create an image. The 
signals from fat and water are at slightly different frequencies and therefore appear 
1 5 shifted relative to each other in an image. 

The fat/water shift is relatively small when a low field, clinical MRI system is 
used. Fortuitously, a much stronger spectroscopic field magnet was used during 
initial efforts at imaging nerve, introducing a much greater displacement of fat in the 
image produced. With the high intensity fat signal shifted away from the nerve, an 
20 enhancement of the nerve's conspicuity was observed. The recognition of this 
enhancement led to the realization that effective neural imaging could, in fact, be 
achieved through the inclusion of fat suppression in system 14. 

Fat suppression apparently enhances the use of conventional MRI systems for 
neurography in several ways. First, the removal of extraneous components reduces 
25 the number of imaged structures to be distinguished. Second, in a fat suppressed 
image a peripheral nerve exhibits a relatively high intensity and will stand out sharply 
against the low intensity space left behind by the suppressed fat. As will be described 
in greater detail below, fat suppression also synergistically increases the apparent 
magnitude of diffusion anisotropy and magnetization transfer effect. 
30 One suitable fat suppression technique involves the use of a chemical shift 

selective (CHESS) pulse sequence, described in detail, for example, in Haase et al. 
NMR Chemical Shift Selective Imaging, 30 PHYS. MED. BIOL. 341-344 (1985). 

As shown in FIGURE 1 1 A, CHESS involves the application of a sequence of 
narrow band rf pulses A, B and C to the excitation coil 62 to selectively excite the 
35 nuclear spins of fat molecules within the region R of the patient being imaged. By 
way of example, three millisecond Gaussian pulses having a minus three dB 
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bandwidth of 600 Hertz may be employed. A sequence of gradient pulses a, b, and c 
is then applied to the three sets of gradient coils 64 and 66 to dephase the excited 
spins, thereby minimizing the contribution of the fat signals to the final image. The 
gradient pulses a, b, and c applied to the orthogonal gradient coil pairs produce, for 
5 example, gradients of five Gauss per centimeter for three milliseconds along the x, y, 
and z-axis, respectively, 

FIGURE 12 illustrates the effect of fat suppression on neurograms produced 
with the MRI system 14. The image provided in FIGURE 12 is of the forearm of a 
rabbit and corresponds to the images of FIGURES 1-5 described above. The darker 

10 portions of the image represent greater image intensity. As shown in FIGURE 12, the 
ulnar nerve UN and median nerve MN are readily identified. 

As an alternative to the use of CHESS for fat suppression, the desired 
suppression may be effected by selective water stimulation. Other suitable 
alternatives include the Dixon technique for fat suppression described in, for example, 

15 Dixon et al.. Simple Proton Spectroscopic Imaging, 153 RADIOLOGY 189-194 
(1984) and also STIR (short tau inversion recovery) described in Improved Fat 
Suppression in STIR MR Imaging: Selecting Inversion Time through Spectral 
Display, 178 RADIOLOGY 885-887 (1991). 

Although in the preferred embodiment fat suppression is combined with other 

20 techniques such as diflfijsional weighting and long T2 processing, fat suppression by 
itself enhances conventional MRI processing sufficiently to generate clinically useful 
neurograms. Similarly, as will be described in greater detail below, other techniques 
employed by system 10 can be used without fat suppression to generate suitable 
neurograms. 

25 ii. Spin-Echo Sequence (Without Diffusional Weighting) 

Having discussed the optional introductory portion of the illustrative pulse 
sequence depicted in FIGURE 11, the n^ phase of the neurography system's 
operation will now be considered. 

In that regard, an rf excitation pulse D, shown in FIGURE 1 1 A, is applied to 

30 coil 62 to tilt the net magnetic moment of the spins by ninety degrees relative to the 
polarizing field, into the transverse plane. The resultant maximum transverse 
magnetization then decays to zero as the spins dephase, A second pulse E, having 
twice the intensity of pulse D, is applied to coil 62 after a delay of one-half the return 
or echo time (TE). This pulse rotates the spins a further 180 degrees and causes a 

35 spin-echo to form as the spins rephase. The spin echo has a maximum amplitude after 
a further delay of TE/2. A spin-echo signal F is, thus, generated in coil 62 at time TE 
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in response to the combined influence of excitation pulse D and refocusing pulse E. 
These steps are depicted in blocks 104, 106, and 108 of FIGURE 9. 

At the same time, the imaging gradients are produced by the orthogonal coil 
pairs 64 and 66 to encode the echo signal F in the usual manner, allowing an MR 
5 image to be constructed, as indicated in block 110. With the sample oriented along 
the z-axis, the "slice select" pulses d, d', and e shown in FIGURE IIC are applied to 
the 2-axis coil pair 64 and 66, to excite and refocus the z-axis slice of interest. The 
"readout gradient" pulses f and f, shown in FIGURE 1 ID, are applied to, for 
exstmple, the x-axis coil pair 64-66 to achieve the desired output that is to be Fourier 

10 '^ ^^x^^^ ^L^'^T\^ "phase encoding" pulses g and g', shown in FIGURE HE, are 
applied to the y-axis coil pair 64 and 66, to control the number of echoes (e.g., 256) 
to be received. The sequence may be used to generate images from contiguous slices 
or regions of the patient. 

As will be appreciated, if the operator indicates (at block 100) that diffiisional 

15 weighting is not required for the generation of a particular image by neurography 
system 10, the pulses shown in FIGURES IIA-IIE define substantially the entire 
spin-echo sequence. Even if diffiisional weighting is to be employed, in the preferred 
embodiment an initial image is generated using only fat suppression for enhancement 
and, as a result, diflfusional weighting is not used during the first performance of the 

20 spin-echo sequence (blocks 104-1 10) for a particular slice. 

Although spin-echo imaging is employed in the preceding embodiment of 
neurography system 10, other techniques can be employed. Suitable alternative 
techniques include, for example, stimulated echo imaging and gradient-recalled echo 
imaging, e.g., echo planar imaging (EPI). Such alternative techniques are described in 

25 Parikh, MAGNETIC RESONANCE IMAGING TECHNIQUES ( 1 992). 

iii. Echo Processing 

In the imaging sequence depicted in FIGURE 1 1 , a series of echo signals F are 
acquired to create a two-dimensional image. For example, at block 112 in 
FIGURE 9, 256 echoes with 256 different phase encoding gradient amplitudes are 

30 used to construct a 256-by-256 pixel image. The data set is then enlarged at 
block 114 by zero filling to produce a 1024-by-1024 matrix of data. As a result, the 
apparent resolution of the final image is increased, making the image clearer. 

Next, the enlarged data set is processed using a 2D Gaussian filter at 
block 116. The filter smoothes the image by attenuating the high frequency 

35 components in the image and, thus, clarifies the delineation of small details without 
altering the relative average pixel intensities over a region of interest. At block 118, 
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the two-dimensional matrix of data then undergoes a two-dimensional Fourier 
transform, which yields an image to be stored. If desired, the image may also be 
displayed on the computer monitor, although in the preferred arrangement this image 
is but one component used in a more extensive analysis performed to generate a 
5 select, enhanced image. 

Once an initial image has been generated, the analysis of the image is initiated, 
as shown in FIGURE 10. At block 122, one or more regions of interest (ROI) within 
the image can be identified. Each ROI may be a single pixel or voxel, or a larger 
region. ROI selection can be performed manually using, for example, a keyboard or 
10 mouse to move a cursor over the ROI on the displayed image. Alternatively, ROI 
selection may be accomplished automatically via a sequential selection of all pixels or 
via an ^^ernal input regarding a particular region from, for example, diagnostic 
system ^ro. 

Next, the average image or pixel intensity within each ROI is computed at 
15 block 124. This average image intensity S can be represented by the following 
expression: 

S = Ao[exp(-TE/T2)][exp(-bD)] (1) 
where Aq is the absolute signal intensity for a particular pixel and the b is the gradient 
factor, determined in accordance with the expression: 

20 b = y2(Gi2)(52)(A.5/3) (2) 

where y is the gyromagnetic ratio, Gj is the polarizing field strength, 5 is the length of 
a diflfiisional weighting gradient pulse, and A is the interval between dif[usional 
weighting gradient pulses. As will be appreciated, in the first iteration before 
diffiisional weighting is employed, the final term of equation (1) is, thus, unity. 

25 To make use of the expressions in equations (1) and (2), the preceding data 

acquisition process is repeated for different values of echo time TE. On the other 
hand, if diffijsional weighting is employed, as described in greater detail below, the 
data acquisition process is repeated for different gradient strengths (controlled by 
adjusting gradient magnitude and/or duration) or gradient orientations. For example, 

30 TEs of 30, 60, 90, and 110 milliseconds, or gradient magnitudes of 0, 3, 5, and 7 
Gauss/centimeter, may be employed. The image intensity S for a particular pixel of 
these multiple images of the same transverse slice for particular values of TE (or b, if 
diffiisional weighting is employed) is available and a linear regression analysis of the 
logarithmic relationship is performed at block 126. 

35 Finally, the value of the apparent T2 relaxation time (or the apparent diffiision 

coefficient D, if diffijsional weighting is employed) is computed for a particular ROI 
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at block 128. These computations provide quantitative assessments of the various 
ROI in the image that are useful in subsequent image processing by other components 
of the medical system 12. 

iv. Gradient Orientation for Diffusion Weighting 
5 In the preferred arrangement, afler the initial fat suppressed image has been 

collected and its ROI characterized, a diffusional weighting analysis is initiated to 
further enhance the neurograms generated by evaluating the diffusional anisotropy 
exhibited by nerve and other tissue. The first aspect of this analysis is the selection of 
the diffusional gradients to be used. 

10 By way of introduction, in one currently preferred embodiment, the analysis 

involves the application of pulsed magnetic field gradients to the polarizing field in 
two or more directions to produce images in which the peripheral nerve is enhanced 
or suppressed, depending upon the "diffusion weighting" resulting from the particular 
pulsed gradient axis chosen. Discrimination of water diffusion anisotropy is then 

15 achieved by subtracting the suppressed image from the enhanced image, in the manner 
described in greater detail below, producing an image depicting only the peripheral 
nerve. 

Most preferably, the magnetic field gradients are applied in mutually 
substantially orthogonal directions. For example, with gradients approximately 

20 perpendicular and parallel to the axis of the peripheral nerve at the particular point 
being imaged, the parallel gradient image can be subtracted from the perpendicular 
gradient image to produce the desired "nerve only" image. 

As will be appreciated, if the axis of the nerve is generally known to the 
operator and its relationship to the referential frame of the MRI system 14 has been 

25 indicated at initialization block 100, the direction of the desired orthogonal diffusional 
weighting gradients can be readily determined. On the other hand, if the axis of the 
peripheral nerve is not known, or if many^nerves having different axes are being 
imaged, the neurography system 10 must employ a system of gradient orientations 
suitable for imaging nerve having substantially any axial alignment. For example, as 

30 will be described in greater detail below, a full three-dimensional vector analysis can 
be used to characterize the diffusion coefficient and provide a nerve image by 
construction based upon a fixed arrangement of diffusion weighting gradients. 

In anatomical regions, such as the upper arm or wrist, it is also possible to 
achieve adequate enhanced isolation of the nerve image by applying only a single 

35 diffusion gradient perpendicular to the axis of the nerve at the site of interest. As a 
result, no subtraction need be carried out to produce the neurogram. The fat 
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suppressed, orthogonally diffusion weighted image can either be processed directly, or 
it can be subject to threshold processing to remove signals of lower intensity 
associated with non-neural tissue, or nerves with different axes and directions of 
travel at the imaging location. 

As will be appreciated, for quicker and more efficient data collection and 
processing, the establishment of diffusion gradients in the polarizing field should be 
responsive to the particular one of the foregoing scenarios that applies to the imaging 
problem at hand. Depending upon the inputs provided at block 100, the system may 
have been advised that (a) only one gradient of known orientation is required, (b) two 
orthogonal gradients of known orientation are required, or (c) two or more gradients 
of unknown orientation are required. 

As indicated in block 130 of FIGURE 10, upon completion of the analysis of 
an image, the system considers whether all of the desired diffusional gradients have 
been applied to the polarizing field during subsequent data acquisition by, for 
example, spin-echo processing, or fast spin-echo processing. Because no diffusional 
gradient was used in the initial fat suppression processing, the answer is initially NO 
and operation proceeds to block 132. 

There, the computer determines whether the operator initially indicated that 
the axis of diffusional anisotropy is known. If the axis is known, a perpendicular 
diffusional gradient is employed, as indicated at block 134. Then, as indicated at 
block 136, a diffusion-weighted spin-echo sequence is performed (modified by the 
inclusion of the diffusional gradient in the manner described in greater detail below) 
and image generated, pursuant to blocks 102-122, before quantification of the image 
data/ occurs at blocks 124-128 to compute D or T2. If the operator indicated at 
initialization that orthogonal diffusion gradients are required for the particular imaging 
problem at hand, this process is then repeated at blocks 138 and 140 for a parallel 

diffusional gradient. - 

If the inquiry performed at block 132 determines that the axis of diffusional 
anisotropy is unknown, operation proceeds to block 142. There an initial difflisional 
gradient is arbitrarily selected, to be followed by a sequence of alternative gradients 
selected for use by the operator when the anisotropic axis is unknown. 

At block 144, using the initial diffusional gradient, a spin-echo sequence is 
performed (modified by the inclusion of the diffusional gradient in the manner 
described in greater detail below) and image generated, pursuant to blocks 102-122, 
before quantification of the imaged data occurs at blocks 124-128. Then, at 
block 146, a test is performed to determine whether the desired number of different 
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diffusional gradients (e.g., three gradients, along the x-, y-, and z-axes) have been 
used. If not, the next diffusional gradient is selected at block 148 and the spin-echo 
sequence, imaging and processing operations are performed, as indicated at 
block 144. This process is then repeated until the desired number of alternative 
5 diffusional gradients have been employed. 

As will be appreciated, additional gradient coils may be provided where 
gradients are desired along axes other than those provided by the locational gradient 
coils. To that end, diffijsional gradient coils may be mounted on a magnetically 
compatible, adjustable track within the bore of the imager to allow gradients to be 

10 repositioned and applied over a substantially continuous range of orientations. 
Similarly, the region to be imaged may be movably supported relative to a fixed set of 
gradient coils to introduce the desired variability in gradient direction. As another 
option, a plurality of different gradient coils may be employed and activated in various 
combinations to effect the desired gradient variations. Alternatively, the results 

15 obtained from a limited number of gradient directions can be processed using a vector 
analysis to estimate the results obtainable with a gradients other than those directly 
available, as described in greater detail below. 

V. Spin-Echo Sequence For Diffusional Weighting 

As noted briefly above, for each of the different diffusional gradients 

20 employed, the spin-echo sequence is repeated, followed by the generation of image 
data and the processing of that data to, for example, quantify the relaxation time T2 or 
diffusion coefficient D. In the preferred arrangement, the use of diffusion gradients 
influences a number of aspects of the spin-echo sequence. 

As shown in FIGURE 1 IF, two pulses h and h', applied to the desired pair of 

25 gradient coils are used to establish a particular diffusional gradient in the polarizing 
field. For an echo time (TE) of 50 milliseconds, the duration (5) of each pulse is, for 
example, 10 milliseconds and their separatioH (A) is 20 milliseconds. In the presence 
of the diffusional gradient, the echo signal F, and therefore the pixel or voxel intensity 
in the image ultimately produced, is made sensitive to the spatial diffusion of water 

30 molecules in the imaged region R. 

In that regard, as indicated above, with the diffusional gradient oriented 
substantially perpendicular to the difflisionally anisotropic nerve, the nerve image is 
enhanced and generally exhibits the highest intensity of various features imaged. This 
phenomena is depicted in FIGURE 13 A, which is an image of the forearm of a rabbit, 

35 corresponding to the diagram provided in FIGURE 1. The ulnar nerve UN and 
median nerve MN are both relatively dark (high intensity) and are easily seen. 
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Altematively, with the diffiisional gradient oriented substantially parallel to the 
diffiisionally anisotropic nerve, the nerve image is suppressed and generally exhibits a 
lower intensity than other features imaged, as illustrated in FIGURE 13B. These 
images can be combined, via a subtraction process described in greater detail below, 
5 to produce an image of the nerve isolated from all other structure. 

To reduce the effect of cross-terms between the imaging gradients and the 
diffusion weighting gradients, the spin-echo sequence illustrated in FIGURE 1 1 is a 
modified version of conventional sequences. More particularly, the readout gradient 
rephasing pulse f shown in FIGURE 1 ID, is placed directly before the acquisition of 

10 echo F in FIGURE 1 1 A, instead of after the slice-selective excitation pulse d, shown 
in FIGURE IIC. However, a consequence of this change was the appearance of 
artifacts in the non-diffiision-weighted images due to an unwanted echo, presumably 
formed from imperfections in the slice-selection pulses d, d*, and e, shown in 
FIGURE lie. To overcome this problem, a second modification of the pulse 

15 sequence was made. Specifically, the phase-encoding gradient was split into two 
sections g and g', and two or four transients (depending upon S/N) were acquired 
with phase cycling. As a result, the remaining cross terms contribute less than three 
percent to the diffusion weighting factor. 

Although fat suppression is not required to take advantage of the image 

20 enhancements available through diffiisional weighting gradients, in the preferred 
arrangement, the fat suppression sequence shown in FIGURES llA and IIB is 
employed prior to the initiation of the diflfiision-weighted spin-echo sequence. As will 
be described in greater detail below, the combination of these techniques generally 
provides an image quality that exceeds that available from either technique 

25 individually. 

The echo F produced using the diffusion weighted pulse echo sequence is 
processed in the manner described above in Connection with blocks 112 through 128 
of FIGURES 9 and 10. With diffusion weighting, the computation of the diffiisional 
coefficient D at block 128 is preferably based upon the analysis of data collected for 

30 different gradient magnitudes. For example, the computation may be based upon 
gradients of 0, 3, 5, and 7 Gauss/centimeter, resuhing in the production of image data 
as represented in FIGURES 14A through 14D, respectively. While fat, bone, 
marrow, skin and vessels are generally absent even at the lower gradients, muscle and 
ligaments drop out at the higher gradients. As suggested previously, the increasingly 

35 stronger gradients may be achieved by increasing gradient duration, rather than 
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magnitude. Alternatively, the iterative data collection process may be performed 
using different gradient directions. 

vi. Image Selection/Production 

Once computer 72 determines, at block 130, that images have been collected 
5 for all of the desired diffiisional gradients, operation proceeds to block 150. If the 
axis of anisotropy is unknov^n, the various diffiisional coefficients D computed for 
each ROI using different gradient orientations are compared at block 150 to identify 
the maximum and minimum values. These coefficients provide a measure, associated 
with each pixel or voxel, of the magnitude of diffiisional anisotropy at that point, 
10 while the anisotropic direction is indicated by the gradient orientation. 

(a^ Subtraction Neurography 

In the preferred arrangement, the images associated with the maximum and 
minimum values of the diffiisional coefficients for a particular ROI are then used in a 
subtraction process, as indicated at block 152. The image associated with the larger 

15 coefficient is produced by a gradient that is more nearly perpendicular to the neural 
axis, enhancing the nerve image, while the image associated with the smaller 
coefficient is produced by a gradient that is more nearly parallel to the axis, selectively 
destroying the nerve signal. When these two penultimate images are then 
mathematically (or photographically or optically) subtracted from one another, a 

20 subtraction neurogram is produced. 

By way of illustration, FIGURE 15 A is an image produced without diffiision 
weighting. FIGURES 15B and 15C then illustrate images produced using parallel and 
perpendicular gradients, respectively. Finally, the subtraction neurogram produced 
when the image of FIGURE 15C is subtracted from that of FIGURE 15B is shown in 

25 FIGURE 15D. 

This "ideal" neurogram is somewhat analogous to a subtraction angiogram (an 
image showing only blood vessels), but sharpSy highlights a nerve rather than a vessel. 
Such an image is particularly useful for confirming the identification of nerves in a 
given imaging plane or space as well as for locating nerve injuries and nerve 

30 compressions. Despite the well known existence of angiograms showing the entire 
vascular pattern in an anatomic region, and despite the existence of MRI techniques 
that could have been applied to the problem of neural imaging techniques, and despite 
the great need for the visualization of nerves, particularly, in isolation, there has not 
previously been any way of creating such neurograms. 

35 Although image subtraction is employed in the preferred arrangement, it is not 

necessary. For example, in some applications of known anisotropy, subtraction is 
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unnecessary and can be foregone in favor of a threshold analysis. Also, the 
subtraction process can be further supplemented, if desired. For example, the output 
of the subtraction process can be divided by the signal information from a fat 
suppressed, T2-weighted spin echo sequence (e.g. using the aforementioned CHESS 
5 technique). 

One potential problem to be addressed by the use of the subtraction process is 
image registration. As will be appreciated, provided that non-neural tissue is 
identically located in both images subjected to the subtraction process, the non-neural 
component will cancel out of the resultant image. On the other hand, if some shift or 

10 other discrepancy in the apparent position of non-neural tissue is introduced into an 
image due, for example, to movement of the subject, cancellation may not occur and 
the nerve may actually be more difficult to identify in the resultant image. 

In one embodiment, acceptable image registration is evaluated prior to 
initiation of the subtraction process. More particularly, the intensities of pixels in one 

15 image are compared to the intensities of corresponding pixels in the second image. 
Pixels of neural tissue are disqualified on the basis of their high diffusional anisotropy, 
assessed via their diffusion coefficients. Unless the intensities of the remaining, non- 
neural pixels fall within a certain range of each other, indicating acceptable image 
registration, subtraction will be inhibited. 

20 fb^ Vector Processing and Three-Dimensional Image Generation 

Up to this point, the output produced is generally in the form of a single two- 
dimensional image, or a series of two-dimensional images that can be related to form 
a three-dimensional image. In a simple form of three-dimensional image generation, 
described in greater detail below, the high S/N ratio of the two-dimensional 

25 neurograms produced by system 14 readily allows the imaged nerve cross-sections to 
be identified and then linked together to form a three-dimensional projection of the 
neural structure. 

As will be appreciated, however, depending upon the neural pattern involved 
and the spatial resolution required, this simplified approach may introduce undesired 

30 discontinuities into the three-dimensional projection, A more sophisticated processing 
scheme employs information about the anisotropic direction of the nerve in each two- 
dimensional image to further enhance the accuracy of three-dimensional image 
projections. The availability of information regarding anisotropic direction is also 
useful in establishing the optimal directions for the gradients used in the difJiisional 

35 weighting analysis described above to produce a two-dimensional image. 
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In that regard, the anisotropic axis of the peripheral nerve being imaged is 
sometimes known to the operator, allowing the operator to input the directional 
information at block 100 and select the best diffiisional gradients for imaging. More 
commonly, however, nerves and CNS neural tracts follow relatively complex paths 
5 and the direction in which the diffusion anisotropy coefficient of the nerve or tract is 
greatest gradually shifts from one plane or axis to another as the nerve or tract curves 
or turns. As a result, one or two arbitrarily oriented, standard gradients may be 
inadequate to provide the desired images. 

Changes in neural direction can be monitored by moving the patient relative to 

10 a fixed set of gradient coils or employing movable diffiisional gradient coils mounted, 
on a track with a non-magnetic drive system, within the bore of the imager to 
adjustably control the orientation of the diffiisional gradients applied to the region of 
interest. By monitoring changes in the ratio of Dpj/Dpr obtained for a given pixel 
using alternative gradient alignments, or for sequential pixels using the same gradient 

15 alignments, changes in neural direction can be estimated and suitable gradient 
directions selected. Alternatively, gradient coils oriented in three planes can be 
simultaneously activated in various combinations to achieve the effect of an infinite 
variety of differently oriented gradients. 

One advantage of attempting to track changes in neural direction is that 

20 parallel and perpendicular gradient information can then be collected and used to 
produce a subtraction neurogram of the type described above. If, however, the 
optimal gradient directions for a given pixel are determined using feedback from 
images generated with repetitively adjusted gradients, processing speed may be 
significantly impaired. 

25 In many cases the well known anatomy of a nerve will permit the use of a 

particular axis orientation in advance. Initial imaging information will provide a 
description of the gross course of the nerve. - A subsequent "informed" approximation 
can optimize the orientation in each slice. This can be useful to insure excellent 
homogeneity of nerve image intensity or to measure the coefficient of anisotropy 

30 along the course of the nerve. 

As a preferred alternative, requiring less mechanical complexity and faster 
processing speed, a technique has been developed for observing diffiisional 
anisotropy, independent of its degree of alignment with any individual gradient axes. 
This process involves the combination of information from anisotropy measurements 

35 obtained along three standard orthogonal axes or using information from multiple 
fixed axes. For example, in the preferred embodiment, a vector analysis is used to 
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produce interpolated images and directional information from the three orthogonal 
diffusion-weighted images described above. 

In that regard, image information is collected from, for example, four "multi- 
slice" sets using a zero diffusion gradient Bq and diffusion gradients B^, By, B^ in the 
5 X-, y-, and z-orthogonal directions, respectively. For each pixel in the image to be 
produced, information concerning the corresponding pixels in the four diffusion 
gradients images is combined to produce a diffusion vector, representative of water 
molecule movement along the nerve fiber in either direction. This vector has a 
magnitude representative of the image intensity of the pixel and a direction 
10 representative of an "effective" diffiision gradient associated with the pixel. 

More particularly, the image intensity Sn of a given pixel in the new image is 
calculated using the following vector equation: 

Sn = vector length = [{S^^ + Sy2 + S^^ySo^]^^^ (3) 
where S^, Sy, and are the image intensities of the corresponding pixels in the 
15 images produced by the B^^, By, and B^ gradients, Sq is the image intensity of the 
corresponding pixel in the image produced by the Bq gradient and is included in 
equation (3) to normalize the resultant image intensity S^. The direction of the 
effective gradient associated with this pixel image includes components 6^^^, G^z, and 
9y2, computed in the following manner: 
20 Gxy = diffusion vector angle between B^ and By = arctan (Sy/S^) (4) 

Gj^ = diffusion vector angle between B^ and B^ = arctan (SJSj) (5) 
t> Gyz = diffusion vector angle between By and B^ = arctan (Sy/S^)!^ ^ O 

The parameters computed in equations (3), (4), (5), and (6) can be used to 
generate images in a variety of different ways. For example, the intensities of the 
25 pixels can be displayed as a "vector length" image. An illustration of a vector length 
CNS image, in which the intensity of the image is proportional to the magnitude of 
is shown in FIGURE 16. ^ 

The image of FIGURE 16 is a brain scan of a monkey (macaca fascicularis) 
weighing 2-2.5 kg, performed using diffusion imaging (spin-echo) on a General 
30 Electric CSI II imager/spectrometer (2 Tesla, equipped with actively shielded 
gradients). The acquisition parameters were: TR = 1000 ms, TE = 80 ms, diffusion 
gradients = 5 Gauss/cm, diffusion gradient duration = 20 ms, diffusion gradient 
separation = 40 ms. Four slices of thickness 4mm were imaged. T2-weighted images 
were used to reproducibly select the diffusion images. 
35 As an alternative to the use of vector length images, arctan images can be 

employed. These images are obtained by establishing the intensity of a pixel in direct 
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proportion to the angular output of one of equations (4), (5), or (6). An illustration 
of an arctan image is provided in FIGURE 17. As shown in this example of a CNS 
neurogram, a select neural tract of interest can be effectively traced and made to stand 
out in isolation from other neural tracts. 
5 When used to evaluate lesions in CNS images of the type shown in 

FIGURES 16 and 17, vector length images will be more sensitive to water diffusion 
changes where all three orthogonal images change in the same way, while the vector 
angle images will be sensitive to changes in anisotropy between two orthogonal 
directions. A CNS lesion caused by experimental allergic encephalomyelitis induced 
10 by myelin basic protein is demonstrated by its departure from the diffiisional 
anisotropy, which appears as vector length decreases and image intensity changes 
accentuated in particular vector angle images (relationship -ttr^IGURES 16 and 17 
-wieiear). 

Alternative forms of vector analysis can also be applied, for example, as 
15 described in Basser et aL, Fiber Orientation Mapping in an Anisotropic Medium with 
NMR Diffusion Spectroscopy, SMRM BOOK OF ABSTRACTS 1221 (1992). 
Similarly, tensor analyses employing tensors of various ranks, as described in Basser 
et al.. Diagonal and Off Diagonal Components of the Self-Diffusion Tensor: Their 
Relation to an Estimation from the NMR Spin-Echo Sigfiai SMRM BOOK OF 
20 ABSTRACTS 1222 (1992), can be used to treat, or transform the coordinates of, MR 
diffiisional anisotropy data. Suitable alternative processing techniques have been 
developed for use in the evaluation of magnetic, thermal, and structural anisotropy 
data. 

Unlike prior art systems, because the non-neural components of the 
25 neurograms produced by system 14 exhibit a relatively low intensity, or indeed 
disappear entirely from the images, the computer 72 is readily able to identify nerve 
locations in the anatomical structure and to Correctly trace the course of the nerves 
between two-dimensional image planes or through a three-dimensional acquisition 
volume. For example, the location of nerves in a given image plane can be detected 
30 by comparing pixel intensity to some threshold level. A three-dimensional image can 
then be formed by linking or projecting the results of these two-dimensional analyses 
over the desired volume. 

As an alternative, the vector information obtained above can be used to track 
continuous serial changes in the direction of maximum anisotropy of a nerve or neural 
35 tract as the nerve or tract travels along its natural course. In that regard, the direction 
of maximum anisotropy for each voxel associated with a nerve is determined and a 
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voxel connection routine, of the type described in Saloner et al., Application of a 
Connected-Voxel Algorithm to MR Angiographic Data, 1 JOURNAL OF 
MAGNETIC RESONANCE IMAGING 423-430 (1991), is then used to link up 
voxels of maximum anisotropy. The resultant plot of the nerve or neural tract 
5 provides enhanced spatial resolution and less discontinuity from one image plane to 
the next. 

As an alternative to the tv^o-dimensional imaging sequences described above, 
it is also possible to carry out the signal acquisition using a "three dimensional" 
imaging sequence of the type described in Frahm et al.. Rapid Three-Dimensional MR 
10 Imaging Using the FLASH Technique, 10 JOURNAL OF COMPUTER ASSISTED 
TOMOGRAPHY 363-368 (1986). The output of this sequence is then processed 
using a three-dimensional Fourier transform to extract the returns from nuclei over the 
volume being imaged. The resultant processing used to compute D for a given voxel 
and to generate, for example, a subtraction angiogram is substantially the same as 

1 5 described above. 

Regardless of the routine employed to project the neural structure in three- 
dimensions, the system 10 may be further programmed to implement the projection by 
referring to known characteristics of the structure. More particularly, once a given 
nerve has been identified in a given two-dimensional image, an "expert" system 10 is 

20 able to predict the occurrence of certain branches and mergers in this structure, albeit 
at unknown locations. This information can then be used to test the plausibility of the 
projection being generated, refining it where necessary. 

vii. Results of Combined Fat Suppression and Diffusion Weighting 

As previously noted, both muscle and nerve exhibit diffusional anisotropy. In 

25 view of the relatively low signal strength of neural components, diffusional analyses 
were not expected to provide clinically useful neurograms. The combined use of fat 
suppression and diffusional weighting has, -^however, been found to be extremely 
effective in providing the desired nerve image enhancement. 

By way of illustration, for a gradient strength of 7 G/cm and an echo time of 

30 50 ms, an nerve image signal intensity (S^) of 17 and a muscle image signal 
intensity (Sm) of 7 were calculated, based upon the difference between signal 
intensities with pulsed gradients oriented perpendicular and parallel to the nerve. A 
nerve-to-muscle contrast parameter R of 2,43 was then computed as the ratio S^/Sm. 
Similarly, a comparison of the apparent diffusion coefficients for diffusional gradients 

35 perpendicular (Dp^) and parallel (Dpi) to nerve and muscle are as follows: 
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Apparent Diffusion Coefficients 00 ^ cm ^sec't 


Muscle Nerve 

Dpr 1-17 0.65 

Dpi 2.18 2.00 

5 Dpl/Dpr 19 3.1 


These results clearly illustrate that the neural components exhibit a far larger relative 
change in intensity than muscle components when subjected to diffusion anisotropy 
analysis. 

An unexpected and apparently synergistic benefit of fat suppression, when 
10 used in combination with diffusional weighting, is that an actual increase in neural 
signal anisotropy is experienced, with the conspicuity of the neural component of the 
image increasing by roughly 250 percent when the fat component is removed. The 
combined increase in nerve conspicuity and reduction in fat interference significantly 
enhances the effectiveness of neural imaging. 
15 Although not entirely understood, there are several potential explanations for 

the synergistic relationship between fat suppression and diffusional weighting. First, it 
appears that fat suppression may increase the apparent diffusional anisotropy of nerve, 
enhancing the utility of diffusional weighting gradients in the detection of neural 
tissue. By way of illustration, 4fl indicated in the following test data^ obtained with the 
20 signal from fat and "short T2" water removed, the intensity of the remaining image 
signal was due largely to anisotropically diffusing water. 

CHESS Applied CHESS Not Applied 




Gradient Direction 

Gradient Direction 

25 

Nerve Imaeed 

ET 

pi Ratio 

EE 

pi Ratio 


Ulnar Nerve 

29 

<8 >3.6 

62 

49 1.3 


Median Nerve 

30 

<8 ' >3.8 

46 

22 2.1 


Muscle 

14 

8 1.8 

18 

12 1.5 


The synergistic role of fat suppression can also be viewed as a demonstration 
30 of a magnetization transfer effect. More particularly, the irradiation of protons on^ 
e.g,. myelin lipids surrounding the nerve by the saturation pulse of the fat suppression 
sequence may allow transfer of the saturation pulse to water molecules in close 
association with the lipid, allowing for very efficient transfer. Subsequently, these 
molecules can exchange into the anisotropically diffusing mobile water pool. 
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vii. Long TE/TR/T t Processing 

As an alternative to the use of diffusional gradients described above, in some 
regions of interest, it is possible to achieve adequate enhanced isolation of the nerve 
image by use of a spin echo fat suppression technique with a relatively long TE (echo 
5 time) or TR (repetition time) to achieve a T2-weighted image. In that regard, after fat 
suppression, the dominant component remaining in the echo F is returned from 
muscle. Because the T2 of peripheral nerve has been measured by the inventors to be 
roughly twice as long as the T2 of muscles, the use of a relatively long TE or TR in 
the spin echo sequence allows the muscular return to be removed. 

10 The basic operation of a neurography system 14 employing this feature 

remains the same as that shown in FIGURES 9 and 10, except that the initialized 
value for TE is extended. In that regard, the operator may be called upon to initially 
consider whether the desired imaging is likely (e.g., neural imaging in a patient's 
limbs) or unlikely (e.g., CNS imaging) to be disrupted by the presence of muscle. If 

15 muscular interference is likely, a relatively long TE of between 50 and 100 
milliseconds or even longer is initialized at block 100, The particular TE or TR 
selected depends upon the degree of T2 weighting desired. Alternatively, the 
system 14 may be programmed to compare the imaging data separately collected 
using long TE processing and diffusional weighting to assess which provides the best 

20 results. 

Illustrations of the results available with long TE imaging are provided in 
FIGURES 18A through 18D, for TE equal to 30, 40, 60, and 100 milliseconds, 
respectively. In the image of the forearm of a rabbit, provided in FIGURE 18D, 
produced with a field strength of 4.7 Tesla, nerves are brighter than any other 
25 structure in the image. The extent of the increased nerve conspicuity is on the order 
of ten-fold, rendering the images clearljp susceptible for use in constructing 
neurograms. As will be appreciated, lesser ©onspicuities on the order of 1 . 1 may also 
be useful. 

The use of extended TE processing had previously been considered unfeasible. 

30 In that regard, as described in Moseley et aL, Anisotropy in Diffusion-Weighted MRI, 
19 MAGNETIC RESONANCE IN MEDICINE 321, 325 (1991), nerve was believed 
to exhibit a relatively short T2 time. Surprisingly, however, measurements have been 
conducted indicating that the T2 of muscle is approximately 27 milliseconds, while the 
T2 of peripheral nerve is approximately 55 milliseconds, providing a factor of two 

35 difference between the two types of tissue. 
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iv Additional Enhancements for N f-iiral Imaging 
{2l\ Vessel Suppression 

In addition to the fat suppression and muscle suppression techniques described 
above, vessel suppression may be employed to improve the neurographic selectivity of 
the images generated by system 14. Due to the brightness of slowly moving blood in 
5 some otherwise useful sequences, vessel suppression has particular value when used 
in connection with long TE sequence neurograms. 

A variety of alternative approaches can be employed to achieve the desired 
blood vessel suppression. For example, in a first embodiment, the blood vessels are 
separately imaged to produce a flow-based MR angiogram, employing phase contrast 

10 or time-of-flight information. The angiogram may be produced using the MRI 
system 14, under separate program instructions, or the auxiliary data collection 
system 22 described below. The angiogram is then subtracted from the neural image 
to provide a neurogram in which blood vessels content is completely suppressed. 

As previously noted in connection with the discussion of image subtraction, a 

15 registration problem occurs when the information to be removed from one image is 
not identically represented with the same intensity and location in the subtrahend 
image. With vessel image information obtained using an angiography pulse-echo 
sequence (or other techniques described below) and neural image information 
obtained using the neural pulse-echo sequence described above, some difference in 

20 vessel intensity in the two images is to be expected. One way of avoiding registration 
error in this situation is to normalize the angiogram to the corresponding neural image 
(i.e., equalize angiogram intensity based upon comparative measurements at a vessel 
pixel identified on the angiogram). 

A second technique used for vessel suppression is to employ a short TE 

25 sequence to produce a first image in which blood vessels are relatively bright and 
nerves are relatively dim. This image is then subtracted from a second image obtained 
using a long TE and exhibiting bright nerves aJid dim vessels. 

A third blood vessel suppression technique involves the administration of an 
intravenous "black blood" contrast agent to the vessels. The agent is preferably (but 

30 not necessarily) of the "blood pool" type including dysprosium-DOTA poly lysine or 
iron oxide type contrast agent. The blood vessels are thereby blacked out by the 
pharmaceutical agent so there is no need for a subtraction step to produce the desired 
vessel-suppressed neurogram. 

Finally, i carefiiUy adjusted water suppression techniques can be used to limit 

35 the contribution of the blood vessels and cerebro-spinal fluid (CSF) to the neural 
image generated by system 14. One such technique is fluid-attenuated, inversion 
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recovery (FLAIR), described in, for example, Bydde et al., Comparison of FLAIR 
Pulse Sequences with Heavily T2 Weighted SE Sequences in MR Imaging of the 
Brain, 185 RADIOLOGY SUPP. 151 (1992). 
fb^ Motion Suppression 

Some of the image processing techniques described above, including, for 
example, the use of diffusion weighting, may be adversely influenced by motion of the 
region being imaged. To limit the introduction of ambiguous or erroneous content 
into the images produced (i.e., motion artifact), several different hardware and 
software features may be employed by the neurography system 14. 


As shown in FIGURE ^J^S; the acquisition of image information can be 
electively carried out with the region of the patient under examination immobilized by 
a splint 156. The splint 156 includes a rigid base 158, made of plastic or some other 
non-ferro-magnetic material. Base 158 is included to provide a fixed frame of 
reference for the region under examination and is designed to be optionally rigidly 
secured within the bore of the imager L As will be described in greater detail below, 
once the neurography system 10 has imaged the region, the reference frame provided 
by base 158 allows another system, like surgical system 28, to operate within that 
reference frame in a known relationship to the imaged nerves. A non-rigid system, 
described below, employs fiduciary markers applied to the skin surface within the 
splint. 

A rigid frame 160, made of plastic or some other non-ferromagnetic material, 
is attached to base 158 and provides structural support for the splint 156. One or 
more fiduciary markers 162, e.g. water-filled beads or linear marker strips, are 
provided on the frame 160 and/or base 158 to allow the relationship of frame 160 and 
base 158 relative to the imaged region to be determined from the images generated. 
In the preferred arrangement, each marker 162 extends the length of splint 156, so 
that it is visible in each cross sectional image? generated. At least one of the markers 
or strips 162 is aligned at a non-zero angle to the x, y, and z axes of the image plane, 
ensuring that its particular location in the image provides a positional reference to the 
splint. 

A sleeve 164, made of a thin film plastic and filled with a conformable 
substance 166, such as water containing gel, silicone, foam, or cobalt-chloride doped 
water, is formed around the frame 160 and includes straps 168 for use in attaching the 
splint to the patient. As shown, the frame 160 and sleeve 164 include a number of 
open regions 170, providing access to select areas of the region under examination by, 
for example, surgical system 28. Alternatively, the base 158 may be used with two 
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sleeves. A first such sleeve provides complete and continuous skin contact for 
imaging, while a second sleeve secures the arm but generally allows access by a 
surgical device. 

A pump 172 is included to allow the fluid to be introduced into sleeve 164 
5 from a reservoir 174 under pressure, forcing the sleeve against the patient's skin and 
immobilizing the region under examination. A release valve 176, allows the fluid 
within the sleeve 164 to return to the reservoir 174, relieving pressure within 
sleeve 164. 

In addition to reducing motion artifacts, the splint 156 performs several other 

10 functions. First, as suggested above, the splint 156 provides a reference frame that 
can be used by other components of medical system 12 to ensure that actions are 
taken in proper relationship to an imaged neural network. Second, the splint 156 may 
be required to keep the region under examination immobilized for the successful 
administration of therapy or performance of surgery by system 12. 

15 A third function of splint 156 is the reduction of edge effects that might 

otherwise be experienced using fat suppression. In that regard, the surface of a region 
under examination (i.e, the patient's skin) presents an abrupt transition in the nature of 
the material being imaged. The fi^^dgion^geneity caused by this tissue-to-air 
interface causes the fat signals in the^^^t^sm?ace adipose tissue to spread out 

20 and/or shift in frequency relative to deeper lying fat surrounding the nerve. The 
desired effect of fat suppression is, however, to suppress the signal from the 
underlying tissue adjacent the neural tissue to be imaged. By employing thin flexible 
polyethylene or other plastic for sleeve 164 and paramagnetically doped water for 
substance 166, these edge effects introduced at the surface of the region under 

25 examination are reduced on the order of 1 .5 to 5 for nerve tissues. 

Additionally, the splint can be specially designed to move a particular body 
region during an imaging series to introduce 4erial stepped repositionings of the limb. 
The movement can be controlled externally by a hydraulic system with fiber optic 
feedbacks to assess repositioning. In this fashion, it is possible to collect a series of 

30 images with the limb in a controlled series of positions. These can be later assembled 
to provide a kinematic view of stress or impingements upon a nerve during motion as 
for the ulnar nerve at the elbow. 

Motion artifact can also be addressed by the software used to control 
neurography system 10. In that regard, to ensure adequate echo amplitude for MRI, 

35 the net magnetic moment generated by one pulse sequence typically must be allowed 
to return near its equilibrium value before the next pulse sequence is initiated. This 

« 
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factor, in combination with the sheer number of sequences typically required for 
imaging, typically causes data collection to occur over a relatively long period of time 
(e.g., on the order of 1 to 20 minutes or longer). As will be appreciated, the 
likelihood that significant patient motion will be experienced during the data collection 
5 process increases in direct proportion to the time required for data collection. 

The software controlling the data collection processed is optimized to reduce 
at least some of the delays contributing to motion artifact sensitivity. In the 
arrangement described above, information from a number of different images may be 
used to selectively produce a final image. For example, the subtraction neurogram is 

10 typically generated on a pixel-by-pixel (or voxel-by-voxel) basis using information 
from two images obtained with orthogonal diffusion gradients. By interleaving the 
image sequence so that data for a given pixel is collected for each diffusional gradient 
before collecting data for any other pixel, the susceptibility of the subtraction process 
to motion artifact is reduced. Similarly, where multiple images are collected at 

15 different gradient strengths to compute the diffusion coefficient D for a given pixel or 
voxel, as part of the gradient selection process, the susceptibility of the computation 
to motion artifact can be reduced by collecting data for all gradient strengths at one 
pixel before data for other pixels is collected. Thus, despite the relatively simple data 
collection process depicted in FIGURES 9 and 10, in the preferred arrangement, data 

20 collection is interleaved by collecting data in several planes at each acquisition rather 
than completely collecting all repetitions for a given plane before proceeding to the 
next. 

Another technique used to provide the desired motion suppression is based 
upon the anticipation of certain periodic sources of motion that can be monitored by, 

25 for example, the auxiliary data collection system 22. For example, depending upon 
the region of the patient under examination, the patients heartbeat and respiration 
may introduce some motion that is not JT suppressed by the splint 156. With 
information regarding the periodicity of these sources available from system 22, 
computer 72 may then adjust the data collection sequence so that the excitation and 

30 echo pulses occur at consistent times relative to the motion introduced by the sources. 

One technique for reducing respiratory motion artifact in MRI is disclosed in 
U.S. Patent No. 4,930,508 (Shimoni et al.). In contrast, the neurography system 10 
can be used with a variety of techniques including mass-spectrometer monitoring of 
carbon dioxide output, fiber-optic observation at chest wall movement, or auditory 

35 monitoring by long tube stethoscope with automated sound analysis. 
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(r ) Fascicle Identification an H Nerve Enhancement 

Another feature of the neurography system 10 is its ability to image individual 
nerve fascicles. For example, when a phased array coil 62, or other high resolution 
MRI system is used with a long TE sequence, individual nerve fascicles appear much 
brighter than the i^nnSftafln^^ineural tissue within the nerve and between the 
fascicles, and the nerve takes on the appearance of a multifascicled structure. 

By way of illustration, neural images depicting fascicles in the nerve of a 
patient having a nerve graft are provided in FIGURES 20, 21, and 22. These images 
were obtained using a 1.5 Tesla MRI system (Signa System 5.2 software release, sold 
by GE Medical Systems) with standard 1 Gauss/cm gradients and a phased array RF 
coil system of the type described above. A "fast" spin-echo sequence with a TR of 
5000 ms, TE of 102 ms, and 8 echo train was used with fat suppression and spatial 
RF pulses for vessel suppression. 

Two axial series of images were produced using a two dimensional Fourier 
transformation. The first series consisted of 24, five mm thick sections, a 512 x 512 
matrix, one mm skip, and one nex (number of excitations). The second series 
consisted of 41, three mm thick axial sections, a 256 x 256 matrix, zero mm skip, and 
two nex. The field of view was 18 cm and acquisition time was 10.6 minutes for both 
series. 

Images from the second series were post-processed by selecting (manually) an 
elliptical region of interest, approximately two cm in diameter, around the sciatic 
nerve in each of the sections. This region of interest was selected to exclude blood 
vessels, without requiring the use of the more analytically complex vessel suppression 

features described above. 

Projectional images were obtained using a maximum intensity projection 
(MIP) algorithm, available as part of System 5.t (IVI) provided by GE Medical 
Systems. The resultant neurograms show the^interface between the tibial component 
of the sciatic nerve and the surgically placed sural nerve graft, with FIGURES/^ and 
^ illustrating the nerve in progressively larger scale and FIGURE 22 illustrating an 
axial projection of the nerve including graft g. As an additional benefit, this imaging 
protocol depresses the signal from tissues within the nerve, between and among the 
fascicles, so that the individual fascicles (f) of the nerve stand out in sharp profile. 

The ability of the neurography system 10 to image fascicles is important for 
several reasons. First, fascicle imaging enhances the diagnostic usefiilness of the 
neurogram because it makes it possible to observe and analyze the internal structure 
of the nerve for evidence of disease. As will be appreciated, this observation and 
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analysis may be performed visually by the operator or automatically, as part of the 
operation of the diagnostic system 24 described in greater detail below. 

Second, this unique internal organization can be used to provide neural 
selectivity and enhancement in the imaging process, even when the conspicuity or 

r 

5 signal intensity of a particular nerve does not permit identification. More particularly, 
blood vessels, lymphatics, lymph nodes and collections of adipose tissue, are often 
similar to nerve in shape, location, and intensity in the cross-sectional images. None 
of these features, however, exhibit the internal fascicular structure of nerve. 

By way of illustration, fascicle identification and nerve confirmation may be 

10 used to distinguish nerve from other structures in an ambiguous image in the 
following manner. First, a thresholding process is used to identify relatively bright 
regions of the image potentially representative of nerve. With the boundaries of these 
regions established, the intensity of the pixels associated with each region is evaluated 
and average image intensities for the regions are computed. 

15 If the intensity of a given pixel within a region is more than some 

predetermined amount below the average intensity, the structure associated with that 
pixel is a potential fascicle, A positive fascicle identification is, however, only made if 
one or more of a plurality of predetermined sequences of such pixel groups 
representative of fascicular structure are identified. For example, a group of at least 

20 3 such pixels must be found which are adjacent each other and bounded on at least 
one side by pixels not satisfying this criteria. 

The results of this analysis can be used to distinguish bright regions associated 
with nerve from those associated with, for example, blood vessels or lymphatics. The 
image intensity of regions not satisfying the fascicular identification parameters may 

25 then be adjusted to zero, effectively eliminating these ambiguous structures from the 
image. Alternatively, a neuroradiologist or other specialist may use this information 
to select a volume of interest which the neifrography system can then render into a 

projection neurogram. 

(d) Miscellaneous Nerve Enhancement Techniques 
30 The various embodiments and features of the neurography system 10 

described above can be modified to incorporate alternative approaches to nerve 
identification and enhancement. 

For example, a magnetization transfer pulse sequence can be employed after 
the fat suppression sequence to enhance neural imaging. Magnetization transfer 
35 involves the excitation of chemically shifted protons with an "off resonance" pulse. 
These protons in a short T2 isotropically diffusing water compartment then exchange 


UOFW\8060AP.DOC 



-39- 


into a long T2 anisotropically diffusing compartment. In doing so, they carry the high 
intensity magnetization signal with them, thus inducing a transfer of magnetization to 
surrounding neural tissue to increase its conspicuity in the image. Nerve may exhibit 
efficient exchange between the off'-resonance, relatively stationary protons in the 
5 myelin sheath and the resonant, mobile protons of axoplasmic water. On the other 
hand, muscle does not exhibit exchange with a large off'-resonant proton pool to a 
comparable degree. The magnetization transfer pulse sequence is designed to exploit 
this differential sensitivity between nerve and muscle by using stimulation methods 
similar to fat suppression to synergistically improve the neurographic selectivity of the 

10 image^n two ways simultaneously. 

Other alternative pulse sequences can also be used. For example, a version of 
steady state free precession (SSFP), as described in Patz et al., The Application of 
Steady-State Free Precession to the Study of Very Slow Fluid Flow, 3 MAG. RES. 
MED. 140-145 (1986), can be used. The SSFP is, however, modified to be included 

15 in an imaging protocol to achieve fat suppression. Similarly, a magnetization 
prepared rapid gradient echo (MP-RAGE) sequence, as described in Mugler et al., 
Three Dimensional Magnetization Prepared Rapid Gradient-Echo Imaging (3D MP 
RAGE), 15 MAG. RES. MED. 152-157 (1990) can be used if modified to improve T2 
contrast. In addition, neural selectivity can be achieved by employing proton fast 

20 exchange rates or T^ relaxation rates. 

Further alternative techniques for generating neurograms employ sequences 
optimized to be sensitive to the slow coherent flow of the endoneurial fluid. These 
sequences provides a unique signal because of the proximal to distal direction of flow 
and because of the slow flow rate which can be monitored by techniques originally 

25 developed to distinguish diffusion from perfusion. Such techniques include, for 
example, velocity compensation by gradient moment nulling as described in Ehman et 
al.. Flow Artifact Reduction in MRI: A Review of the Roles of Gradient Moment 
Nulling and Spatial Pre-saturation, 14 MAG. RES. MED. 293-307 (1990) and 
Moran, A Flow Velocity Zeumatographic hiterface for NMR Imaging, 1 MAG. RES. 

30 IM. 197-203 (1982). 

3. Medical System Construction and Operation 
As noted previously, the neurography system 10 is one component of a 
broader medical system 12. The remaining components of system 12 are described in 
greater detail in the following sections. These components provide information to, 
35 and process information from, neurography system 10 in accordance with software 
instructions executed by, for example, a host processing system 32 or the processing 
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systems of individual components of the system 12 to achieve a variety of functions 
beyond the imaging of peripheral nerve. 

a. Auxiliary Data Collection Svstem 
The auxiliary data collection system 22 may take any one of a variety of 

5 different forms. For example, as suggested above, system 22 may be designed to 
collect supplemental information regarding structure present within the images 
produced by system 10. Examples of such systems include a secondary MRI system, 
employing conventional pulse-echo sequences suitable for use in angiography or STIR 
sequences to show areas of high muscle signal due to denervations or fianctional loss 

10 of the muscle; an X-ray imaging system suitable for use in generating image data of 
bone and/or tissue, a PET scanning system for showing the progress of an axionally 
transported pharmaceutical a^^t; or a CT system for collecting contrast agent 
lymphography data. The splint^is also formed with fiduciary markers visible using CT 
and MRI (e.g., iodine contrast material in water), allowing the information from 

15 systems 10 and 22 to be integrated. 

The supplemental information may be used to suppress structural content in 
the image and provide greater neural selectivity. For example, an angiogram may be 
used to remove vessel image content from the neurogram or to distinguish nerve and 
vessels on the basis of color. Alternatively, because non-neural structure is generally 

20 absent in the image anyway, the additional information may be employed to add 
specific structures, such as blood vessels, back into the neurogram unambiguously. 
This process allows alternative structure to be readily differentiated using different 
colors to display information from different sources. As will be appreciated, the 
addition of structure into a MR image viewed to assess neural structure was virtually 

25 unthinkable with prior art systems due to the low signal content of neural return 
components. 

An alternative type of data collecfion system 22 is employed to collect 
information about the patient for use in controlling the operation of the neurography 
system 10, rather than modifying its output. Examples of such systems include 
30 conventional heart rate and respiration monitors, used to time the data collection 
sequencing of system 10 relative to the heart rate and respiration of the patient. 

A final type of data collection system 22 of interest is one designed to collect 
supplemental information about the neural network. For example, system 22 may be 
constructed to produce an output indicative of nerve conduction velocity (NCV), 
35 including the approximate location of a change in NCV or the NCV response to 
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magnetic stimulation. Information from evoked potential electrodes or magnetic 
SQUID detectors might also be collected and integrated for a multi-input display. 


The diagnostic system 24 is selected to process the image neurograms and 
other information (such as D and T2) provided by neurography system 10 to provide 
an attending physician with, for example, diagnoses of neural anomalies. 
Alternatively, system 24 may assist the physician in making a diagnosis, (^assessing 
the need for, or likely success of, surgery. In one embodiment, system ^may be 
employed simply to confirm or question the physician's diagnoses. 

By way of illustration, one region in which problematic neural disorders 


commonly occur but are difficult to ^mgfto^s-Ts the spinal canal. As shown in the 
cross sectional view of one vertebra provided in FIGURE 23, the region of interest 
exhibits a relatively high physiological complexity. The illustrated structures include a 
herniated disc (HD), compressed left spinal root (LSR), spinous^gro^s (SP), anulus 
fibrosus (AF), nucleus pulposus (NP), autonomic ganglion (AG),/yentral root (LVR), 
ventral ramus (VRA), transverse process (TP), dorsal ramus (DRA), dorsal root 
ganglion (DRG), facet (F), dorsal root (DRO), extradural fat (EF), root in cauda 
equina (RCE), dural sac (DS), and cerebrospinal fluid (CF). 

In this diagram, the two features of primary interest are the left spinal root 
(LSR) and the left ventral root (LVR), which are both in risk of compression from the 
herniating disc (HD). Both nerves are traveling through extradural fat (EF) but are 
surrounded by bone, which could impair observation by an X-ray based technique. 
Both nerves are also near the strong water signal of the disc (HD), cerebrospinal fluid 
(CF) in the dural sac (DS), and other inflamed tissue (which often diminishes image 
resolution and quality in generally used magnetic resonance techniques). 

The left spinal root (LSR) and left ventral root (LVR) of diagnostic interest 
are small relative to the numerous large anatomic structures nearby. Also, these roots 
are nearly perpendicular to each other. This common imaging problem can be 
addressed by the use of a neurography system 10 programmed to employ fat 
suppression, followed by pulsed diflaision gradients oriented to enhance either the left 
spinal root (LSR) or left ventral root (LVR), so that each can be clearly seen in a 
given image. If either root is compressed, its image will demonstrate physical 
distortion or an imprint due to the compression, which may manifest itself as a change 
in structure or signal intensity between the two sides of a compression. 

As will be appreciated, due to the selectivity and resolution of the neurograms 
produced by system 10, they can be evaluated by a physician to diagnose any neural 


b. Diagnostic System 
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abnormalities present. In addition, the image produced by system 10 can be analyzed 
by diagnostic system 24 to detect, for example, evidence of compression or 
inflammation and provide the appropriate diagnosis. 

The operation of system 24 depends, in part, upon the condition to be 
5 evaluated. In one embodiment, the operator initially views a two- or three- 
dimensional image generated on a cathode-ray tube (CRT) display, included with 
system 24 and uses a cursor to identify a particular imaged nerve to be evaluated. 
The operator may also input the particular types of anomalies to be detected. 

The system 24 then determines the boundaries of the imaged nerve in each of 

10 the two dimensional images available, using a thresholding process. These boundaries 
can then be compared from one image to the next to look for discontinuities or 
changes in shape associated with a particular anomaly of interest. For example, if this 
analysis were used with a nerve severed in an accident, the nerve might disappear 
entirely from certain images in which it would otherwise be expected to appear. The 

15 system 24 is able to readily identify such regions and provide the physician with 
precise locational information regarding the anomaly. 

Similarly, the physician may be interested in conditions associated with less 
pronounced changes in nerve boundary or intensity. The system 24 is readily able to 
provide outputs indicative of the average intensity of a bounded neural area, as well as 

20 the size and shape of the bounded area on an image-by-image basis. This information 
can then be used to detect anomalies such as compressions. 

In one arrangement, a cursor can be used to initialize a reference boundary of 
interest on the CRT (e.g., associated with a "normal" neural cross section) for use by 
system 24. System 24 then compares the actual boundary of the nerve in subsequent 

25 images to the reference boundary to locate and quantify the extent of neural 
compressions. This quantification of an anomaly then allows the physician to monitor 
the recovery of the nerve and assess the effecliveness of any therapy being provided. 

Another approach that may be employed by system 24, is based upon the 
apparent increase in T2 exhibited by injured nerve. More particularly, an initial "long 

30 T2" analysis or diffusion weighted image can be performed to image all neural 
structures. Then T2 can be extended to roughly 100 milliseconds to image only those 

nerves that are injured. 

An additional approach for use in the manual evaluation of, for example, bone 
fractures and injuries to joints, involves the analysis of an image in which the fat 
35 component is selectively demonstrated and remaining tissues suppressed. This 
approach emphasizes the appearance of skin, adipose collections, and of bone (in 
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many locations) due to the presence of marrow. When such an image is collected and 
assigned a color other than that used to display nerve, the two images can be shown 
transparently in the same three dimensional construction. As a result, the physician is 
provided with useful information regarding the physical relationship between nerves 
5 and bones. This information is most important in the evaluation and treatment of 
bone fractures and joint injuries. 

In some applications, a contrast agent may be used to synergistically highlight 
the anomaly of interest. Alternatively, because nerves appear bright and isolated in an 
image, it may be more informative to selectively black out one of the nerves by means 
10 of administering an intraneural pharmaceutical contrast agent. 

In addition to analyzing the output of neurography system 10, the diagnostic 
system 24 may also provide feedback to system 1 0 to control the pulse sequences 
used and the type of information produced. For example, where sites of nerve 
compression, section, laceration, or fibrosis are imaged, the alteration in endoneurial 
15 fluid flow and in axoplasmic flow are readily detected by monitoring the increase in 
signal intensity when T2-based, or other, neurographic sequences are used. 

Although not described in detail herein, a variety of different diagnostic 

applications are contemplated including: 

1. The demonstration of a patient's peripheral, cranial, and autonomic 

20 nerve and nerve plexus anatomy. 

2. The demonstration of a patient's spinal root anatomy, particularly the 
cervical, thoracic and lumbar spinal roots and nerves where they pass through fat at 
the foramina through which they exit the spinal canal. 

3. The demonstration of a patient's spinal root anatomy within the lumbar 
25 canal where the roots pass through quantities of extradural fat. 

4. The examination of a patient's cranial nerves for compressions by 
vessels or other structures which could cause trigeminal neuralgia (Vth nerve), 
hemifacial spasm or Bell's palsy (Vllth nerve), essential hypertension (Xth nerve) or 

other cranial nerve syndromes. 
30 5. The demonstration of nerve, plexus or root compressions or injuries in 

a patient, where abnormal changes in the direction, position, or other diflRisional 
properties are caused by an injurious process, such as nerve transection, demyelinating 
diseases, neuritis, multiple sclerosis, peripheral neuropathies and crush injuries, as well 
as the monitoring of the regrowth of nerves. 


UOFW\8060AP.DOC 


-44- 


6. The determination of the location of tumors or other masses within the 
spinal cord where it is useful to know the position of cortico-spinal motor tracts or 
other functional white matter long tracts relative to some abnormality. 

7. Demonstrating the anatomy of the optic nerve, an extension of the 
5 brain, where it passes through the peri-orbital fat^spther fat on its route to the retina. 

8. Tract tracing within the brain to provide useful images for study by 
radiologists, surgeons or physicians and, in particular, for identification of the location 
of areas of 'eloquent cortex' such as the motor strip, or speech-related areas. This 
method involves the spatial identification of relevant areas of the thalamus or internal 

10 capsule and then following projecting tracts to the area of interest on the cortical 
surface, or identifying regions of interest by reference to their connections with other 
areas on the cortical surface. For example, speech cortex projection tracts can be 
followed from areas known to be involved in speech production to (and through) 
other areas where an injury or stroke may have blocked proper nerve function. 

15 9, Tracing of nerves as they pass through tumors of low diffusional 

anisotropy, such as the passage of the Vllth nerve through an acoustic neuroma to 
permit a surgeon to know the location of the nerve in or near the tumor and so to 
have the ability to avoid the nerve during surgery on the tumor. 

10. Application of diffusion anisotropy imaging for the evaluation of 
20 diffuse axonal injury, as may occur in head injury. 

11. The evaluation of bone fractures and joint dislocations or 
dislocation/fractures in which surgical planning, management and fixation would 
benefit from knowing the course of the nerve in the region of the abnormal anatomy. 

c. Therapeutic Svstem 

25 As noted previously, a therapeutic system 26 is also employed to process 

information from neurography system 10 or other components of medical system 12 
to better effect the administration of therapy-."to the patient. For example, system 26 
may be a drug-delivery system or a current-stimulation system that employs feedback 
from neurography system 10 to regulate its operation. In this fashion, more precise 

30 nerve conduction velocity (NCV) or evoked potential tests can be done using 
neurographic data to place stimulating or recording electrodes. For therapy, tract 
information could aid in the placement of transplant tissue or for lesions of areas of 
abnormal activity that might cause tremor in the thalamus. 

d. Surgical Svstem 

35 The surgical system 28 employs neurographic information from system 10 to 

influence any one of a variety of surgical operations that may be performed. The 
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information obtained may be used to avoid neural paths during surgery or to confirm 
the location and nature of neural surgery required. The operation of surgical 
system 28 may be automatically controlled in response to feedback from system 10 or 
may manually controlled by a surgeon based upon his or her review of the information 
5 provided- 

In one embodiment, the region of the patient that is to undergo surgery is 
placed in the splint 156 described above. The open regions 170 of splint 156 need to 
be designed and positioned to ensure that splint 156 will not interfere with the surgical 
system 28 during the performance of an operation. With the splint 156 applied, image 

10 data is then collected via the neurography system 10. As noted previously, the 
processing system 16 of neurography system 10 provides numerical coordinates, in 
three dimensions, describing the position of the nerves along their courses with 
reference to the splint base 158 and fiduciary markers 162. 

Depending upon the nature of the operation to be performed, outputs from the 

15 auxiliary data collection system 22 may also be required. For example, if system 28 is 
employed to operate on bone within the region imaged, system 22 may be called upon 
to generate a fat selective image of the bone, or the patient may be brought to a C-T 
scanner for preparation of a bone image. The MRI splint 156 will be worn while this 
additional information is collected, but additional markers (e.g., chalk or iodine 

20 solution for CT X-rays) are required to extract locational information from the 
secondary image and, hence, provide the requisite registration between the two 
images generated. 

The image information is loaded into the memory of a surgical system 
processor 178, shown in FIGURE 24. As will be described in greater detail below, in 

25 the preferred arrangement, processor 178 is programmed to guide surgical operations 
performed in a coordinate system that is referenced to the image coordinate system. 
The base 158 of splint 156 is secured to a pl^form 180 included in surgical system 28 
to provide a fixed relationship between the coordinate systems used in the image and 
by system 28. The coordinate systems are then linked by processor 178 using a 

30 computer model of three-dimensional space. Confirmatory X-rays may be taken 
conveniently during the procedure to assure correct positioning of the markers. 

An articulated surgical arm 182 is coupled to the platform 180 and has a 
stylus 184 (e.g., a surgical apparatus, such as a focused laser beam or a drill) provided 
at its free end. The arm 182 can be moved electrically or pneumatically to any select 

35 point, or along any select path, defined relative to the operating environment in 
response to outputs from processor 178. The position of the arm 182 can, thus, be 
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tracked via the control outputs applied by the processor 178. As will be appreciated, 
a separate coordinate-based or laser-based positioning system may be employed if 
desired. 

In the preferred embodiment, the position of the stylus relative to the imaged 
5 neural and other networks is illustrated on a system display 186 during the course of a 
surgical operation. The surgeon may manually guide the stylus 184 during the 
operation via, for example, a joystick, electronic glove, or other input device 188, 
visually monitoring the position of the stylus relative to anatomical structure. This 
visual feedback may be based simply upon a comparison of the known positional 

10 relationship of the stylus to the previously collected image. 

Alternatively, it is possible to obtain visual confirmation using imaging 
feedback data collected in real time. For instance, with the use of high-speed MRI 
data collection sequences, such as echo planar imaging described in Worthington et 
al.. The Clinical Applications of Echo Planar Imaging in Neuroradiology, 32 

15 NEURORADIOLOGY 367-370 (1990), it is possible to rapidly update images. 
When the resulting images are displayed, the surgeon may observe the progression of 
an appropriately labeled, non-magnetic probe into the body in real time. If a slower 
image collection process is employed, the probe or device is advanced in steps as a 
series of images are taken. In either case, the neurographic image provides the 

20 surgeon with apparent vision of sensitive neural tissue inside opaque, solid body 
structures, in much the same manner as fluoroscopy, but while also providing 
information regarding neural paths. 

As an alternative to requiring the surgeon to control the operation of 
system 28 during surgery, a computer-guided, stereotaxic, or fiduciary system may be 

25 employed. In that regard, the surgeon may provide processor 178 with input 
identifying the nature of the operation to be performed, including the stylus path and 
operation appropriate for the surgery to be performed. These steps can be performed 
with the arm 1 82 disengaged, allowing the surgeon to simulate the operation and view 
the stylus path on the image, before the surgical procedure is confirmed. Once 

30 confirmed, the processor 178 can then be instructed to automatically guide the 
arm 182 over the desired path during the actual surgical operation. 

The use of surgical system 28 has a number of important advantages over the 
"neurally blind" surgical methods currently employed. For example, because nerves 
are readily imaged, the surgeon is better able to assess any neural conditions that 

35 might require treatment or alter the surgical plan. In addition, because the position of 
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the stylus 182 relative to nerve is readily imaged and can be confirmed before an 
operation is performed, accidental intrusion of the stylus upon neural paths is avoided. 

Although a splint 156 is employed in the embodiment described above to 
provide a link between the referential frames of the neurogram and surgical system 28, 
5 it is not mandatory. For example, particularly in regions that are relatively 
unsusceptible to motion artifact, fiduciary markers can be applied directly to the body 
(e.g., on the head or face when nerves of facial sensation or movement are involved, 
or adjacent the iliac crests and lumbar vertebral spinous processes when lumbar nerve 
roots are involved). 

10 The use of a computer-guided surgical system 28 employing such fiduciary 

markers is believed to be of particular importance in cervical, thoracic or lumbar spine 
surgery. In that regard, system 28 will eliminate the problem of doing a "good" 
operation but at the wrong level, e.g., inadvertently decompressing the lumbar 3/4 
root when the symptomatic compression to be relieved was actually at the lumbar 4/5 

15 root. For spine work, the original image can be collected with a strip of fiduciary 
markers taped to the patient's back and independently marked in a manner that can be 
sensed by system 28 to allow location of the stylus 182 during surgery. 

As will be appreciated, the various components of the surgical system 28 can 
be altered in a variety of manners. For example, the stylus 182 may include a surface 

20 detector of electrical fields or a magnetic detector of nerve activity, constructed to 
detect the activity of nerves. Examples of such devices include a somatosensory 
evoked potential or magneto-encephalography system. As a result, the detection of 
nerves offered by embodiment of the stylus 182 allows the position of the stylus 
relative to nerve determined by reference to the neurogram to be confirmed. 

25 One application for surgical system 28 that is of particular importance is in 

surgery of the neck. Surgery of this type includes, for example, a carotid 
endarterectomy to remove stroke producing^plaque from the internal carotid artery, 
an anterior cervical discectomy to relieve a cervical root compression, or an operation 
for cancer in the neck. One of the most common complications of such surgery is the 

30 accidental crushing or transection of the recurrent laryngeal nerve, possibly resulting 
in permanent paralysis of one or both sets of vocal cords. Optimally, a preoperative 
neurographic image is used to demonstrate the course of the recurrent laryngeal 
nerve, allowing the surgeon to more effectively avoid it or at least identify and protect 
it during surgery. 

35 Neurographic guidance can be also used for percutaneous needle biopsy of 

lesions, or for the placement of more elaborate percutaneous systems such as 
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ultrasonic or other mechanical devices used to remove tissues. By way of illustration, 
such operations include discectomies, the introduction of laser/suction systems, the 
placement of RF lesioning devices used in procedures such as gangliolysis of the fifth 
cranial nerve, the placement of probes to carry out deep tissue localized drug 
5 administration, diathermy, cryotherapy, or other physical or mechanical techniques. 
Neurographic guidance may also be used to control the passage of rigid endoscopes 
through solid tissues or to guide the placement of directable flexible endoscopes. 

Yet another important application of surgical system 28 is in the use of CNS 
neurograms to guide stereotactic surgery in the brain. Currently, tissue structures 
10 visible by virtue of their or T2 MRI are used to guide stereotactic surgery. In 
contrast, CNS neurograms provide information concerning the connections or relation 
of specific tracts of interest, which may travel in or among other tracts from which 
they cannot be differentiated by means of conventional tissue-based images. 

e. Training and Development Svstem 
15 The training and development system 30 may take any one of a variety of 

forms designed to process information collected from the neurography system 10. In 
one embodiment, neurographic images are collected fi-om a plurality of patients to 
produce an anatomical atlas of normal and abnormal neural paths for reference by 
surgeons and others. Images obtained from a patient can be compared to the 
20 clinically known population in the atlas to rapidly identify anomalous nerve courses in 
a patient set to undergo surgery for some condition. As a result, the surgeon may be 
able to modify his or her technique to reduce the risk of injury to nerves which happen 
to be in the field of surgery. Similarly, a neurographic map of an individual patient's 
skin and cutaneous nerves can be used to help the surgeon plan incisions that avoid 
25 the very common complication of accidental transection of cutaneous nerves in the 
course of routine surgical incision to reach structures below the skin. 

In another embodiment, the training and development system 30 may be 
designed to assess the effectiveness of the programming employed by neurography 
system 10 and may provide feedback to system 10 to regulate its operation and 
30 enhance the quality of the neurograms generated. More particularly, once a sequence 
able to positively identify nerve has been employed, alternative sequences can be 
employed and their results compared to the confirmed method. As a result, a 
collection of techniques can be established and programmed into neurography system, 
along with the conditions under which each sequence offers the best performance. 
35 Another alternative training and development system 30 may be employed to 

assess the effectiveness of intraneural, pharmaceutical contrast agents designed to 
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help in the diagnosis of nerve compressions. More particularly, such a training 
system 30 is used to quantify the image contrast produced using different contrast 
agents to image known neural anomalies. As a result, system 30 is able to identify 
those agents providing the best results for particular neural imaging problems. 
5 Yet another embodiment of training and development system 30 allows 

information from neurograms to be used in the design of any of a variety of products. 
For example, the neurograms produced by system 10 provide information that is of 
great advantage to designers of ergonomic furniture, high gravity air or space craft 
seats, specialized body suits, boots, and various kinds of electronic or electric medical 

10 equipment, which can be best used when the positions of neives can be precisely 
located in advance. The system 30 incorporates information regarding neural paths 
from system 10 into the computational processes involved in designing such 
equipment to provide' enhanced product performance. 

As one illustration, in the ergonomic design of a chair, system 30 would be 

15 programmed to ensure that the primary support provided by the chair does not 
coincide with a neural path in the relevant customer population. This can be done by 
superimposing the neural network of a sitting person upon a mathematical model of 
the chair, identifying the primary points of support, and generating flags on the display 
for any support points that are within some predetermined distance of a nerve. As a 

20 result, the chair design can be manipulated to avoid neural compressions. 

In another application, the system 30 can be used to control electr^iic _ 
(y prosthesis. More particularly, the information from system 10 can be used to focatedT^ 

electronic detectors adjacent, for example, a severed nerve to detect neural activity 
associated with the limb replaced by the prosthesis. The detected activity of the nerve 

25 is then used to control the prothesis. 

4. Non-neural Imaging Applications 
In principle, selective imaging of any^other object or subject may be effected 
using the MRI techniques described above, if that subject or object exhibits 
characteristics corresponding to the neural characteristics described above. For 
30 example, objects exhibiting diffusion anisotropy in any part thereof can be imaged 
using diffusional weighting. Thus, in medicine, for example, the cardiovascular/ 
^ system could also be imaged in this fashion and the technique can also be employed to 

examine, for example, rock strata and plants, if they exhibit diffusion anisotropy. 

5. Conclusion 

35 The lack of a suitable method for creating a distinct image of the nerves has 

been a great hindrance to physicians, surgeons, athletic trainers, and pain treatment 
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specialists. Although, previously, it has sometimes been possible to make a nerve 
stand out from immediately surrounding structures, the unique ability of system 10 to 
make the nerve stand out from all other structures represents a significant advance. 
The sensitivity of the system 1 0 allows even the smallest nerves to be accurately 
5 identified and linked to form three-dimensional projections of a neural network. The 
neurographic information can be collected rapidly, without requiring contrast agents. 

While the preferred embodiment of the invention has been illustrated and 
described, it will be appreciated that various changes can be made therein without 
departing from the spirit and scope of the invention. 
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